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control  techniques  for  these  types  of  systems.  An  Hx  controller  is  synthesized  for  a 

complex  analytical  model  with  997  inputs,  936  outputs,  and  332  states.  To  accomplish 
this,  a  new  technique  for  model  reduction  using  Zernike  polynomials  was  developed. 
The  Hx  controller  was  able  to  achieve  a  minimum  15  Hz  control  bandwidth.  The 
previous  integral  controller  was  unable  to  meet  the  10  Hz  bandwidth  requirement.  The 
Hx  design  process  used  was  validated  on  a  simpler  adaptive  optics  testbed.  The 
experimental  verification  also  showed  that  the  robust  control  techniques  outperformed  the 
classical  control  techniques  in  the  presence  of  disturbances. 

The  significant  contributions  are  a  Zernike  polynomial  method  for  model 
reduction,  robust  controller  synthesis  for  a  complex  adaptive  optics  analytical  model,  and 
experimental  verification  on  an  AO  testbed.  Although  the  robust  control  design  is  more 
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I.  INTRODUCTION 


A.  MOTIVATION 

Optical  systems  have  long  been  used  for  imagery  and  intelligence  gathering  for 
both  civilian  and  military  purposes.  Since  the  first  recorded  use  of  telescopes  in  the  early 
17th  century,  the  goal  has  been  to  improve  the  technology  that  would  give  the  user  the 
advantage  of  increased  resolution  at  greater  distance.  Some  of  the  early  pioneers  in  the 
field  such  as  Galileo  and  Newton  experimented  with  refractive  (lens)  and  reflective 
(mirror)  systems.  As  the  technologies  allowing  remote  optical  sensing  evolved,  military 
commanders  on  the  battlefield,  as  well  as  naval  commanders  at  sea,  found  new 
applications  for  them.  Naturally,  the  military  application  of  these  tactics  evolved  with  the 
advent  of  air  and  space  vehicles  from  which  to  conduct  reconnaissance.  The  trend  has 
continually  moved  towards  smaller  angular  resolution  at  increased  distance. 

As  shown  in  Equation  (1.1),  the  diffraction  limited  angular  resolution  using 
Raleigh  criterion,  as  denoted  by  aR,  is  a  function  of  the  wavelength  of  the 

electromagnetic  radiation,  approximately  600  nanometers  for  visible  light,  and  the 
diameter,  D,  of  the  light  gathering  aperture 

ccR=1.22A/D  (1.1) 

The  angular  diffraction  limited  resolution  is  usually  measured  in  radians, 
arcseconds,  or  some  equivalent  measure.  Since  the  wavelength  being  observed  is  usually 
fixed,  the  aperture  diameter  is  the  only  remaining  variable  factor.  Obviously,  there  are 
physical  size  limitations  to  which  the  primary  optic  in  a  telescope  can  be  manufactured. 

The  most  obvious  limitation  arises  from  the  physical  process  of  machining  the 
lens  or  mirror.  Ideally,  a  reflective  mirror  would  have  a  parabolic  shape  to  minimize  the 
amount  of  aberration  resulting  from  the  reflection  process.  In  practice,  however, 
spherical  mirrors  are  much  easier  to  machine,  giving  rise  to  one  of  the  most  common 
optical  aberrations,  appropriately  known  as  spherical.  Additional  constraints  on  the 
diameter  of  the  primary  optic  arise  in  the  support  structure  used  to  encase  or  hold  the 
optic  stationary.  Large  telescope  optics  will  often  be  quite  rigid,  a  design  characteristic  to 
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prevent  sagging  or  warping  of  the  optical  surface  due  to  gravitational  or  other  inertial 
forces.  Mechanisms  to  support  and  move  these  heavy  optics  can  themselves  become 
quite  large  and  cumbersome. 

This  support  structure  problem  has  less  consequence  for  ground-based  observing 
platforms,  but  for  airborne  platfonns,  weight  becomes  more  important.  This  realization 
becomes  immensely  important  when  considering  space-based  platforms  where  the  cost 
per  pound  to  orbit  is  prohibitively  high.  Space  platforms  are  perhaps  even  more  limited 
by  size  constraints  than  weight  constraints.  The  basic  shape  of  every  launch  vehicle  that 
launches  spacecraft  into  orbit  is  cylindrical.  Even  the  space  shuttle  orbiter,  whose 
exterior  is  not  cylindrical,  has  a  cylindrical  payload  bay  that  carries  the  spacecraft  aloft. 
The  diameter  of  the  rocket  body  becomes  a  physical  limitation  for  the  size  of  the  primary 
optic  of  a  space-based  telescope.  Considering  the  distance  involved  in  orbital  altitudes, 
this  diameter  is  certainly  less  than  that  needed  to  achieve  the  angular  resolution  desired 
by  image  analysts. 

Despite  these  challenges,  the  attraction  of  space-based  imagery,  and  in  particular, 
high  altitude  imagery,  is  quite  appealing.  Beyond  the  immediate  advantages  of  a  wider 
field  of  view  and  higher  vantage  point,  a  high  altitude  orbit  also  results  in  lower  orbital 
speeds,  which  translates  to  increased  dwell  time  over  a  given  part  of  the  earth.  This  can 
be  taken  to  the  extreme  of  the  geostationary  orbit,  which  would  maintain  continual 
presence  over  a  point  above  the  equator  and  have  nearly  a  third  of  the  surface  of  the 
planet  in  its  field  of  view  at  any  time.  The  tradeoff  comes  in  the  distance  to  geostationary 
altitude.  Whereas  a  low  earth  orbiting  satellite  might  orbit  a  few  hundred  kilometers 
above  the  surface,  geostationary  spacecraft  require  an  altitude  of  approximately  36,000 
km.  The  effect  of  this  distance  is  illustrated  in  Figure  1  and  Equation  (1.2),  where  the 
angular  resolution  is  equal  to  that  from  Equation  (1.1),  and  R  represents,  in  this  case, 
distance  from  the  orbital  altitude  to  the  surface  of  the  earth. 
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RaR  =  GD  (1.2) 

where  GD  is  the  ground  resolution. 

Although  increasing  altitude  does  not  decrease  the  angular  resolution  per  se,  it 
does  increase  the  minimum  separation  needed  to  differentiate  two  objects  on  the  surface. 
Said  another  way,  space  telescopes  at  geostationary  altitudes  require  much  larger 
diameter  primary  optics  to  differentiate  the  same  level  of  ground  detail  as  low  earth 
orbiting  satellites.  Clearly,  a  way  is  needed  to  combine  the  benefits  of  better  image 
quality  obtained  from  lower  orbits  with  that  of  improved  dwell  time  from  higher  orbits. 
Multiple  solutions  exist  and  are  not  limited  to  individual  spacecraft.  A  concept  of 
operations  could  be  developed  which  utilizes  a  constellation  of  spacecraft  at  varying 
altitudes  to  maximize  resolution  and  dwell  time  over  particular  areas  of  interest.  Another 
solution  involves  achieving  primary  optic  diameters  larger  than  would  fit  inside  a 
conventional  payload  fairing,  by  segmenting,  and  deploying  on  orbit.  These 
segmentation  techniques  will  be  the  topic  of  this  dissertation. 

Traditional  telescopes  utilize  a  spherical,  continuous  face  sheet  mirror.  Its  size  is 
constrained,  however,  by  the  aforementioned  limitations.  If  instead  of  a  continuous  face 
sheet,  the  primary  optic  can  be  made  of  a  series  of  smaller  segments,  these  segments  can 
be  folded  or  packaged  in  some  way  that  the  segments  are  deployed  or  assembled  into  a 
large  diameter  mirror  after  the  rocket’s  payload  fairings  have  been  jettisoned  and  the 
spacecraft  is  on-orbit.  The  concept  is  not  unlike  that  which  is  done  with  spacecraft  solar 
arrays,  which  are  often  folded  like  an  accordion  before  launch  and  unfolded  once  in 
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space.  To  be  clear,  though,  the  on-orbit  assembly  of  mirror  segments  is  vastly  more 
complex  than  that  of  solar  arrays,  communications  antennas,  or  any  other  spacecraft 
structures  deployed  on-orbit. 

The  added  complexity  stems  partly  from  the  fact  that  mirror  segments  are 
extremely  delicate  components  with  sensitive  optical  coatings,  but  primarily  from  the 
precise  alignment  criteria  required  of  the  segments.  Recall  that  the  wavelength  of  visible 
light  is  on  the  order  of  600  mn.  With  wavelengths  this  small,  even  a  slight  misalignment 
of  mirror  segments  can  cause  aberrations  to  degrade  image  quality  beyond  what  is  useful. 
This  problem  is  so  daunting  in  fact,  that  a  separate  area  of  research  is  being  devoted 
entirely  to  the  problem  of  initial  alignment  for  segmented  optics  [1].  Once  the  coarse 
initial  alignment  has  been  successfully  accomplished,  only  fine  adjustments  will  be 
necessary  to  account  for  remaining  optical  aberrations.  More  on  those  aberrations  will  be 
discussed  later.  This  research  will  focus  only  on  the  finer  surface  corrections  after  initial 
alignment  has  been  completed. 

Numerous  folding  and  deploying  scenarios  could  be  considered  for  segmented 
optics.  One  important  design  criterion  is  the  electromagnetic  radiation  wavelength 
collected  through  the  aperture.  The  emphasis  here  will  be  on  visible  wavelengths.  A 
current  example  of  a  space  telescope  under  construction  that  will  utilize  segmented  optics 
is  NASA’s  James  Webb  Space  Telescope,  an  infrared  astronomy  telescope  planned  for 
launch  in  2013.  An  artist’s  concept  of  this  spacecraft  is  shown  in  Figure  2.  Notice  the 
gold-colored  primary  optic,  which  is  comprised  of  18  hexagonal  segments.  The  launch 
configuration  of  this  spacecraft  calls  for  the  three  leftmost  and  three  rightmost  segments 
to  be  folded  back  in  order  to  fit  within  the  payload  fairings.  Two  hinged  motions  will 
deploy  those  six  segments  along  with  a  host  of  other  deployments  from  the  launch 
configuration  involving,  but  not  limited  to,  the  struts  supporting  the  secondary  optic  and 
the  blue  and  red  sun  shield.  Although  the  primary  optic  deployment  may  seem  rather 
simple,  it  is  important  to  remember  that  the  remaining  twelve  segments  are  still 
individual  segments,  and  vibrations  and  stresses  of  the  launch  environment  will  make  for 
a  very  challenging  initial  alignment  problem. 
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Figure  2  James  Webb  Space  Telescope  From  [2]. 


Another  challenging  aspect  of  the  space  telescope  observation  platform  involves 
the  structural  properties  of  the  primary  optic.  Spacecraft  designers  are  constantly 
attempting  to  decrease  the  mass  of  the  onboard  components  to  reduce  the  costs  of  placing 
them  into  orbit.  Large,  rigid,  monolithic  optics  work  fine  for  most  terrestrial 
applications,  but  their  weight  becomes  prohibitive  for  spacecraft  operations.  Whether 
segmented  or  a  continuous  face  sheet,  space  telescope  optics  are  most  economical  when 
lightweight.  Being  lightweight  causes  particular  challenges  in  the  control  and  alignment 
of  mirror  segments  or  other  means  of  mirror  shape  control.  The  natural  frequencies  will 
be  lower,  resulting  in  control  and  structure  interactions. 

Both  segmented  and  continuous  face  sheet  mirrors  will  have  actuators  attempting 
to  correct  the  mirror  shape  due  to  measured  aberrations.  If  the  bandwidth  controlling 
these  actuators  overlaps,  or  even  falls  close  to  the  natural  frequencies  of  the  structure, 
then  the  controller  may  actually  exacerbate  the  aberrations.  Typically,  large  ground- 
based  mirrors  have  their  lowest  natural  frequency  well  over  1000  Hz  so  control 
bandwidths  below  100  Hz  have  no  chance  of  exciting  natural  frequencies.  With 
segmented  mirrors,  such  as  the  one  researched  in  this  dissertation,  the  lowest  natural 
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frequency  can  be  around  30  Hz,  which  becomes  a  concern  for  the  control  system 
designer.  Additionally,  for  space  platforms,  the  natural  frequencies  are  often  tightly 
packed  and  have  very  little  damping.  The  robust  control  techniques  researched 
demonstrate  an  innovative  way  to  design  controllers  for  these  lightweight  segmented 
space  telescopes. 

B.  ADAPTIVE  OPTICS  SYSTEMS 

The  term  adaptive  optics  refers  to  any  optical  system  in  which  some  optical 
component;  be  it  a  mirror,  lens,  or  some  other  device,  undergoes  adaptation  or 
modification  to  improve  the  resulting  wavefront.  The  preceding  discussion  on 
controlling  the  shape  of  the  primary  optic  on  space  telescopes  fits  the  general  description 
of  an  adaptive  optics  (AO)  system.  To  be  clear  on  the  nomenclature,  the  adaptive  optic  is 
that  component  of  the  system  such  as  the  mirror,  which  is  altered;  whereas  adaptive 
optics  usually  refers  to  the  entire  system. 

The  optical  wavefront  is  a  perpendicular  cross  section  or  slice  of  a  propagating 
beam  of  photons,  which  is  measured  by  a  wavefront  sensing  device  such  as  a  camera. 
For  example,  a  beam  of  coherent  light  from  a  laser  propagates  along  a  cylindrical  path. 
At  that  point,  taking  a  sensor  measurement  at  any  location  or  at  any  instant  of  time  yields 
a  representation  of  the  wavefront.  Ideally,  the  wavefront  from  a  high-quality  laser  would 
be  planar,  or  flat.  In  reality,  aberrations  cause  the  wavefront  to  be  non-planar  and  wavy. 

Aberrations  in  optical  systems  come  from  a  variety  of  sources.  Light  traveling 
through  the  vacuum  of  space  would  be  the  most  un- aberrated  light  possible.  As  the  light 
propagates  though  a  dynamic  atmosphere,  or  any  other  medium,  it  becomes  aberrated  to 
an  observer  on  the  other  side  of  that  medium.  Other  aberrations  arise  from  imperfections 
in  optical  components.  Even  the  human  eye  is  not  immune  to  aberrations,  hence  the  need 
for  corrective  lenses  or  eye  surgery.  For  machined  optics,  the  imperfections  result  from 
uneven  coatings,  imperfect  shapes,  or  impurities  in  glass  among  others.  Light  passing 
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through  even  one  optical  component  has  become  aberrated  beyond  its  original  form.  A 
moving  or  vibrating  platform  can  also  lead  to  wavefront  aberrations,  which  is  of 
particular  concern  for  spacecraft. 

AO  systems  come  in  a  variety  of  styles  and  designs,  but  there  are  usually  three 
components  common  to  all.  Those  components  are  a  wavefront  sensor,  an  adaptive 
optical  component,  and  a  control  computer.  This  is  not  unlike  a  standard  control  system 
that  consists  of  a  sensor,  control  algorithm,  and  actuators  influencing  a  plant.  Each  of 
these  portions  will  be  discussed  subsequently.  A  generalized  schematic  is  shown  in 
Figure  3.  This  research  focuses  exclusively  on  the  design  of  new  and  innovative  control 
algorithms.  The  optical  components  and  wavefront  sensors  used  on  the  experimental 
testbed  are  commercially  available,  off-the-shelf  pieces  of  hardware.  Although  the 
emphasis  here  is  on  controls,  adaptive  optics  is  truly  a  multi-disciplinary  field  and  draws 
on  expertise  in  physics,  structures,  controls  engineering,  and  electrical  engineering  to 
name  a  few. 


It 


Figure  3  Generalized  AO  System  From  [3]. 
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1. 


Wavefront  Sensors 


A  wavefront  sensor  takes  an  optical  measurement  of  some  object  of  interest  and 
represents  that  information  as  a  set  of  discrete  data  points.  One  of  the  most  common 
wavefront  sensors  for  AO  systems  is  the  Shack-Hartmann  wavefront  sensor.  It  consists 
of  two  basic  parts,  a  camera,  and  a  lenslet  array.  The  camera  operates  like  any  CCD  or 
CMOS  camera.  Incident  light  excites  pixels  in  an  array  and  the  light  intensity  is 
converted  to  voltages,  which  is  output  along  with  the  location  of  pixels  that  were  excited. 
The  other  main  element  of  the  Shack-Hartmann  wavefront  sensor  is  the  lenslet  array  at 
the  front  entrance  pupil  of  the  camera.  To  understand  its  purpose,  first  envision  a  card 
that  is  placed  at  the  entrance  pupil,  blocking  all  light  from  reaching  the  camera  sensor. 
Now  poke  a  series  of  holes  in  the  card  in  a  geometric  pattern,  such  as  a  square  matrix. 
Light  from  the  incident  wavefront  will  pass  through  these  holes,  undergo  diffraction,  and 
strike  the  CCD  sensor.  Now  replace  the  holes  with  an  equal  number  of  small  lenses 
(known  as  lenslets)  which  abut  up  against  each  other  such  that  there  is  no  remaining 
surface  of  the  card.  The  entire  blockage  now  consists  of  a  grid  of  lenslets. 

Each  of  these  lenslets  is  identical  and  will  focus  the  incoming  light  on  the  CCD 
behind  it.  The  distance  from  the  lenslet  array  to  the  CCD  is  the  focal  length  of  the 
lenslets.  For  a  planar  incoming  wavefront,  the  location  of  the  focused  spots  on  the  CCD 
will  be  evenly  spaced  in  the  same  geometric  pattern  as  the  lenslet  array. 

If,  however,  the  incoming  wavefront  is  not  planar,  then  the  resulting  spots  on  the 
CCD  will  not  be  evenly  spaced.  Whereas  a  planar  wavefront  will  come  to  focus  directly 
behind  each  lenslet,  an  aberrated  wavefront  will  cause  the  spot  to  come  to  focus  at  some  x 
and  y  displacement  on  the  CCD.  Figure  4  shows  an  example  of  an  x  displacement  from 
an  aberrated  wavefront;  a  y  displacement  would  be  nonnal  to  the  surface  of  the  page. 
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Figure  4  Lenslet  Array  Spot  Formation  From  [3]. 

By  measuring  the  actual  location  of  the  spots  vs.  the  expected  location  for  a  planar 
wavefront,  an  estimate  of  the  true  shape  of  the  incoming  wavefront  can  be  obtained.  This 
information  is  then  used  by  the  control  computer  to  detennine  what  action  must  be  taken 
to  correct  the  aberration. 

The  Shack-Hartmann  is  only  one  example  of  a  wavefront  sensor  for  an  AO 
system.  A  phase  diversity  sensor  is  another  common  wavefront  sensor  that  attempts  to 
determine  wavefront  shape  based  on  the  difference  in  image  intensity  from  typically  two 
cameras  placed  at  different  locations  in  the  beam  path  [4],  Both  the  analytical  model  and 
experimental  testbed  described  in  this  research  use  Shack-Hartmann  wavefront  sensors. 

2.  Adaptive  Optical  Components 

Once  the  actual  shape  of  the  wavefront  has  been  detennined,  the  control 
algorithms  are  needed  to  detennine  how  best  to  make  the  wavefront  as  planar  as  possible. 
First,  however,  it  is  best  to  understand  the  actual  methods  available  to  apply  a  correction 
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to  the  wavefront.  This  is  accomplished  through  use  of  an  adaptive  optic.  Perhaps  the 
simplest  type  of  adaptive  optic  to  start  with  would  be  a  Fast  Steering  Mirror  (FSM).  The 
reflecting  surface  of  an  FSM  is  a  normal  flat  mirror.  The  difference  is  in  the  mounting 
mechanism  that  supports  the  mirror  in  the  system.  Imagine  a  pair  of  orthogonal  axes 
along  the  surface  of  the  mirror.  These  axes  are  commonly  referred  to  as  tip  and  tilt  axes. 
By  making  small  rotations  about  these  axes,  the  beam  of  light  the  mirror  is  reflecting  can 
be  steered.  This  ability  to  alter  the  orientation  of  the  mirror,  and  therefore  the  beam  path, 
allows  an  FSM  to  be  classified  as  an  adaptive  optic  and  an  optical  system  containing  an 
FSM  is  likewise  classified  an  AO  system.  An  FSM,  however,  corrects  only  first  order 
wave  aberrations. 

One  common  source  of  such  aberrations  in  an  optical  beam  is  jitter  caused  by 
structural  vibrations.  Since  these  vibrations  are  somewhat  oscillatory  in  nature,  if  the 
wavefront  sensor  is  fast  enough  to  identify  their  frequency,  an  FSM  should  be  able  to 
apply  an  appropriate  phase  difference  to  null  out  the  jitter  aberrations  in  the  optical  beam. 
In  this  example,  the  system’s  actuators  affect  a  change  in  the  FSM’s  tip/tilt  orientation. 

For  higher  order  wave  aberrations,  a  more  complex  adaptive  optic  such  as  a 
deformable  mirror  is  used.  Unlike  the  FSM,  which  has  a  continuous,  flat,  face  sheet 
surface,  the  surface  of  a  deformable  mirror  can  be  adjusted  to  take  on  other  desired 
shapes  that  are  not  flat.  This  is  accomplished  with  a  set  of  actuators  on  the  back  surface 
of  the  mirror.  Piezoelectric  actuators  are  one  type  and  can  be  bonded  to  the  back  surface 
such  that  the  actuator  can  be  lengthened  or  shortened  to  “push”  and  “pull”  the  surface  of 
the  mirror  into  the  desired  shape.  Another  type  consists  of  electrostatic  actuators  where  a 
voltage  potential  creates  a  deformation  of  the  mirror  surface  instead  of  the  motion  of  a 
physically  attached  rod. 

The  purpose  behind  deformable  mirrors  is  to  create  a  mirror  shape  that  exactly 

matches  the  inverse  of  the  measured  wavefront.  Therefore,  if  the  wavefront  is 

determined  to  have  a  particular  shape,  then  by  inverting  that  shape  and  applying  it  to  the 

mirror,  the  wavefront  can  be  made  planar.  The  success  of  this  approach  will  depend  on 

several  factors  including  having  enough  lenslets  in  the  wavefront  sensor  to  accurately 

characterize  the  wavefront,  a  control  algorithm  capable  of  quickly  determining  a  solution 
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before  the  wavefront  changes,  and  enough  actuators  to  create  the  inverse  wavefront  with 
high  enough  fidelity.  Technology  has  increased  to  the  point  where  a  large  number  of 
lenslets  and  actuators  can  be  packed  into  a  relatively  small  space  on  their  respective 
components.  This  is  advantageous  in  that  it  allows  higher  fidelity  shapes  to  be  measured 
and  created.  The  disadvantage  lies  in  the  next  topic  of  discussion,  the  control  laws 
governing  the  process. 

These  two  brief  descriptions  of  adaptive  optical  components  only  scratch  the 
surface  of  the  level  of  complexity  that  can  exist  in  AO  systems.  Remember  that  for  high 
altitude  space  based  telescopes,  the  optics  will  most  likely  need  to  be  segmented  to 
provide  the  degree  of  resolution  required.  A  telescope  comprised  of  even  just  two 
segments  presents  a  significant  challenge  for  the  design  engineer.  Future  telescopes  will 
be  comprised  of  several  individual  segments  requiring  precise  alignment,  where  each 
segment  is  its  own  deformable  mirror  requiring  face  sheet  control.  The  analytical  model 
studied  in  this  research  is  just  such  a  system. 

3.  Control  Algorithms 

Adding  lenslets  and  actuators  to  an  AO  system  not  only  increases  complexity  but 
also  lends  itself  to  coupling  between  actuators.  Whereas  historically,  control  laws  treated 
AO  systems  as  a  static  problem,  and  to  date,  adaptive  optics  has  been  a  field  of  study 
largely  belonging  to  the  field  of  astronomy  where  AO  systems  are  used  to  correct  for 
ground  based  observatories  peering  at  distant  starlight  through  a  turbulent  atmosphere. 
Using  either  a  naturally  occurring  or  artificially  created  guide  star,  the  control  computer 
can  determine  the  atmospheric  aberration  in  the  guide  star’s  wavefront  and  use  that 
information  to  compensate  for  the  atmosphere  with  the  optics.  This  can  be  done  in  either 
a  real-time  or  in  a  post-processing  manner.  For  real-time  correction,  the  computer  has  to 
compute  and  apply  the  corrections  before  the  atmosphere  above  the  telescope  changes. 
This  process  is  made  much  simpler  by  the  telescope  mirrors  being  large  rigid  structures 
where  the  control  bandwidths  are  well  below  the  natural  frequencies,  so  the  structural 
dynamics  can  be  ignored  in  the  control  design.  For  lightweight  space  based  telescopes, 
besides  having  to  contend  with  natural  frequency  concerns,  there  is  increased  probability 
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that  a  large  piston  motion  by  a  particular  actuator  will  have  an  unintended  effect  on  the 
mirror  surface  in  the  neighboring  actuator’s  space.  This  effect  is  magnified  in  light  of  the 
fact  that  on  current  deformable  mirrors,  large  numbers  of  actuators  are  packed  into 
relatively  small  areas.  If  simply  ignored,  this  coupling  will  have  an  adverse  effect  on 
system  perfonnance.  The  structural  flexibility  also  cannot  be  ignored  in  the  control 
design. 

Since  the  design  of  robust  control  algorithms  for  AO  systems  is  the  purpose  of 
this  research,  this  brief  introduction  for  AO  control  algorithms  will  be  sufficient  for  now. 
For  a  more  thorough  treatment  of  AO  systems,  the  reader  is  referenced  to  [5],  considered 
a  primary  resource  in  the  field.  Specifics  of  the  exact  challenges  faced  and  the  solutions 
to  those  challenges  will  be  addressed  throughout  the  remainder  of  this  dissertation. 

C.  OBJECTIVES 

1.  Analytical 

The  analytical  research  has  the  following  objectives: 

1.  Apply  robust  control  techniques  for  fine  surface  control  of  a  complex 
flexible  space  mirror  to  minimize  wavefront  aberrations. 

2.  Develop  improved  model  reduction  techniques  with  the  emphasis  of 
minimizing  wavefront  aberration. 

3.  Achieve  the  desirable  control  bandwidth. 

2.  Experimental 

The  experimental  research  has  the  following  objectives: 

1 .  Apply  robust  control  to  an  adaptive  optics  testbed. 

2.  Compare  the  performance  of  robust  control  and  currently  used  classical 
control  techniques. 
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D.  CURRENT  STATE  OF  RESEARCH  IN  THE  FIELD 

Although  the  majority  of  AO  system  controllers  have  been  of  a  classical  nature, 
there  has  been  some  initial  research  into  advanced  robust  control  techniques  for  AO 
systems.  To  date,  the  author  is  not  aware  of  any  operational  systems  that  employ  a  robust 
controller  based  on  Hx  techniques.  The  work  of  Frazier  et  al.  has  looked  at 

incorporating  some  of  these  techniques  onto  testbed  systems  [6],  [7],  [8],  [9].  Their  work 
has  been  done  on  simple  testbed  setups  quite  similar  to  the  Spacecraft  Research  and 
Design  Center  adaptive  optics  testbed  at  the  Naval  Postgraduate  School  used  in  this 
research.  The  adaptive  optic  in  their  research  was  also  a  37  channel  defonnable  mirror. 

Kun  Li  et  al.  have  investigated  Hx  control  techniques  for  segmented  telescopes 
[10],  but  their  approach  involves  using  centralized,  decentralized,  and  overlapping 
control  architectures  to  handle  the  complexity  of  large  numbers  of  states,  inputs,  and 
outputs.  The  technique  presented  here  will  instead  focus  on  model  reduction  techniques 
to  synthesize  the  Hx  controller.  Other  works  have  begun  to  apply  Hm  techniques  to 

large  ground  based  telescopes  that  have  more  flexibility  than  previous  generation 
telescopes,  similar  to  what  space  based  telescopes  will  experience  [11],  [12],  [13]. 

The  general  theory  of  Hx  control  techniques  is  much  more  developed  than  the 

application  to  AO  systems.  Several  papers  and  textbooks  have  been  written  that  have 
advanced  these  techniques  on  benchmark  problems  [14],  [15],  [16],  [17].  These  works 
have  not  been  specific  to  segmented  space  telescope  applications,  however.  This 
research  will  adapt  these  techniques  for  complex  space  telescope  models. 

The  most  thorough  work  on  Hx  control  techniques  for  lightweight  segmented 

optics  has  been  done  by  Carrier  [18].  His  research  looked  at  modeling  and  robust  control 
techniques  for  the  Advanced  Structure/Control  Integrated  Experiment  (ASCIE)  testbed,  a 
segmented  optical  system  also  used  by  [10].  In  physical  appearance,  the  ASCIE  is  quite 
similar  to  the  space  telescope  model  used  in  this  research.  It  contains  hexagonal  mirror 
segments  for  the  primary  optic  and  a  lightweight  truss  support  structure.  Unlike  the 
model  used  for  this  dissertation,  however,  the  ASCIE  is  not  an  adaptive  optics  system.  It 
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contains  only  24  sensor  inputs  versus  the  936  sensor  measurements  used  here,  and  only 
18  actuators  versus  the  997  used  here.  Also  because  of  some  of  the  computational 
limitations  at  the  time,  the  control  design  was  broken  into  subsystems.  Processor 
capabilities  have  matured  enough  to  negate  the  need  to  form  separate  subsystems; 
however,  a  model  reduction  is  still  required.  Nonetheless,  the  modeling  and  advanced 
control  techniques  developed  by  Carrier  form  the  foundation  for  much  of  the  application 
in  this  dissertation. 

Specifically,  the  work  on  the  ASCIE  demonstrated  that  stochastic  controller 
design  methods  were  impractical  due  to  sensitivities  to  unmodeled  dynamics,  modal 
frequency  uncertainty,  and  modal  damping  uncertainty.  Worst  case,  or  Hv  ,  methods 

yielded  a  10%  stability  margin  to  modal  frequency  uncertainties.  Furthermore,  the 
research  developed  the  technique  of  using  input  and  output  multiplicative  uncertainties 
for  loopshaping  purposes. 

This  research  will  take  the  foundational  work  from  these  sources  and  adapt  it  for 
use  on  a  significantly  more  complex  segmented  space  telescope  model.  Specifically,  the 
techniques  will  be  adapted  for  AO  systems  and  the  requirement  to  perfonn  model 
reduction  prior  to  controller  synthesis.  The  robust  controller  synthesis  process  used  will 
be  experimentally  validated  on  an  AO  testbed.  A  testbed  for  the  analytical  model  was  not 
available,  however,  so  the  control  techniques  will  be  demonstrated  on  a  simpler  setup  and 
suite  of  equipment. 

E.  DISSERTATION  OUTLINE 

Chapter  II  of  this  dissertation  will  give  detailed  descriptions  of  the  complex  space 
telescope  model  along  with  the  AO  testbed  and  its  components.  Chapter  III  will  give  an 
overview  of  robust  control  and  the  factors  considered  in  the  design  of  robust  controllers. 
Chapter  IV  will  describe  the  model  reduction  techniques  necessary  for  the  synthesis  of 
the  robust  controller.  Chapter  V  will  apply  the  model  reduction  and  robust  controller 
synthesis  theory  to  the  analytical  model  and  show  the  perfonnance  when  subjected  to  an 
input  disturbance.  Chapter  VI  will  describe  the  same  process  for  the  experimental  results 
on  the  AO  testbed.  Final  analysis  and  conclusions  will  be  given  in  Chapter  VII. 
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F.  CONTRIBUTIONS 

The  present  research  has  contributed  to  the  state-of-the-art  in  the  following  areas: 

1.  An  Hx  robust  control  technique  has  been  used  for  surface  control  of  a 
flexible  space  telescope  model  including  structural  dynamics. 

2.  Robust  control  techniques  have  been  able  to  provide  a  desired  control 
bandwidth  that  was  not  possible  with  classical  control. 

3.  A  new  model  reduction  technique  using  Zemike  polynomials  and  singular 
value  decomposition  has  been  developed.  This  technique  has  been  found  to 
provide  better  performance  to  reduce  the  wavefront  aberration  in  the  space 
telescope  model. 

4.  A  robust  control  technique  has  been  applied  experimentally  to  the  adaptive 
optics  testbed  of  the  Naval  Postgraduate  School  Spacecraft  Research  and 
Design  Center.  The  robust  control  technique  has  been  found  to  be  superior  to 
classical  control  in  response  time  and  disturbance  rejection. 
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II.  TESTBED  AND  MODEL  OVERVIEW 


A.  SEGMENTED  SPACE  TELESCOPE  MODEL 

1.  Model  Purpose 

The  analytical  model  researched  in  this  dissertation  is  a  state  space  model  for  a 
space  telescope  with  a  segmented  primary  optic.  It  was  developed  and  provided  by  an 
industry  partner  as  part  of  a  government  program  to  test  viability  and  technology 
readiness  for  future  space-based  telescopes.  The  model  represents  an  actual  system  built 
for  test  purposes,  but  not  as  a  flight  test  model.  It  involves  the  latest  state-of-the-art 
technologies  for  space  telescope  systems.  The  analytical  research  conducted  in  this 
dissertation  does  not  involve  the  actual  hardware,  only  the  simulation  model. 

The  model  is  restricted  to  the  optical  payload,  and  does  not  include  the  supporting 
spacecraft  bus.  At  a  basic  level  it  is  not  unlike  other  space  telescopes  with  a  Cassegrain 
configuration  in  that  it  has  a  primary  mirror  that  reflects  the  incoming  light  back  to  a 
secondary  mirror,  which  in  turn  reflects  the  light  back  through  a  hole  in  the  center  of  the 
primary  mirror.  A  suite  of  science  and  calibration  sensors  is  located  behind  the  primary 
optic.  The  James  Webb  Space  Telescope  shown  in  Figure  2  also  has  this  basic 
configuration. 

Due  to  the  nature  of  government  contracts  and  funding,  this  program  was 
somewhat  constrained  in  scope.  For  example,  the  control  system  used  standard  classical 
control  methods  with  notch  filters.  The  nature  of  the  system,  however,  lends  itself  to 
research  into  the  feasibility  of  robust  controllers.  This  is  just  one  of  several  examples  of 
additional  research  that  could  be  accomplished  on  this  model. 

2.  Model  Complexity 

The  segmented  space  telescope  model  is  extremely  complex.  It  contains  332 
states,  997  inputs,  and  936  outputs.  This  is  the  main  reason  that  model  reduction  is 
necessary  before  a  robust  controller  can  be  synthesized.  The  332  states  come  from  166 
modes  identified  in  a  Finite  Element  Analysis.  The  inputs  correspond  to  actuators  for 
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mirror  surface  control  and  the  outputs  are  sensor  measurements.  Most,  but  not  all,  of  the 
sensor  measurements  are  of  the  wavefront.  By  means  of  comparison,  the  segmented 
space  telescope  model  investigated  in  [10],  which  is  actually  quite  similar  in  appearance 
and  function  to  the  model  investigated  in  this  research,  has  70  states,  18  inputs,  and  18 
outputs.  The  Giant  Segmented  Mirror  Telescope  in  [12]  has  200  states  and  91  inputs. 
The  AO  testbed  in  [7]  has  37  inputs  and  162  outputs.  The  benchmark  problems  in  [14] 
and  [15]  typically  have  on  the  order  of  10  states.  Therefore,  as  can  be  seen,  the  scope  and 
complexity  of  this  model  is  much  greater  than  that  of  previous  studies. 

For  this  model,  the  primary  mirror  is  comprised  of  six  identical  hexagonal 
segments.  They  are  arranged  in  a  circular  pattern  such  that  in  the  center  is  another 
similarly  sized  hexagon,  only  vacant,  allowing  the  light  reflected  off  the  secondary  optic 
to  pass  through  to  the  science  instruments.  The  positioning  of  the  individual  segments  is 
accomplished  by  six  coarse  control  actuators  and  three  fine  control  actuators  per  segment, 
for  a  total  of  54  positioning  actuator  inputs.  Each  segment  also  has  a  grid  of  156  face 
sheet  actuators,  for  a  total  of  936  shape  control  inputs,  making  each  segment  a 
deformable  mirror.  Besides  the  adaptive  primary  optic,  this  model  also  has  another 
adaptive  optic,  a  Fast  Steering  Mirror  (FSM),  located  in  the  vicinity  of  the  science 
instruments,  and  is  used  to  remove  global  tip  and  tilt  jitter  motions.  There  are  an 
additional  seven  inputs  related  to  the  operation  of  the  FSM  thus  accounting  for  the  997 
total  model  inputs. 

The  model  outputs  similarly  come  from  several  different  sources.  The  majority 
come  from  a  Shack-Hartmann  wavefront  sensor.  Each  of  the  six  segments  has  6 1  Shack- 
Hartmann  lenslets  associated  with  it.  Each  of  these  lenslets  provides  two  output 
measurements,  an  x  slope,  and  a  y  slope.  This  results  in  732  total  Shack-Hartmann 
outputs.  Figure  5  shows  a  diagram  of  all  six  hexagonal  segments  and  the  individual 
lenslets  associated  with  each  segment.  The  lenslets  are  numbered  for  one  of  the 
segments,  which  will  be  important  for  later  discussion.  The  other  segments  are  numbered 
sequentially  in  a  counter-clockwise  wise  direction  and  the  individual  lenslets  in  each 
segment  are  numbered  similarly  if  each  segment  is  rotated  to  the  top  position. 
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Figure  5  Mirror  Segment  and  Lenslet  Orientation. 


The  model  also  has  a  phase  diversity  wavefront  sensor  providing  1  piston,  1  tip, 
and  1  tilt  measurement  per  segment  for  18  total,  and  a  jitter  sensor  providing  2  total 
tip/tilt  measurements.  The  remaining  sensors  on  the  model  are  18  gap  sensors,  three  per 
segment,  which  measure  the  distance  separating  adjacent  segments.  In  sum,  there  are  770 
total  sensor  measurements.  The  remaining  166  outputs  are  the  position  states.  Table  1 
provides  a  summary  of  the  model  inputs  and  outputs. 
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Inputs 


Sou  rce 

Number 

per  segment 

Face  Sheet  Actuators 

936 

156 

Fine  Control  Actuators 

18 

3 

Coarse  Control  Actuators 

36 

6 

Fast  Steering  Mirror  Actuators 

7 

Total 

997 

Outputs 


Sou  rce 

Number 

per  segment 

Shack-Flartmann  Slopes 

732 

122 

Phase  Diversity  Sensor 

18 

3 

Gap  Sensor 

18 

3 

Jitter  Sensor 

2 

Position  States 

166 

Total 

932 

Table  1  Segmented  Space  Telescope  Input/Output  Summary. 

B.  ADAPTIVE  OPTICS  TESTBED 
1.  Layout 

The  Spacecraft  Research  and  Design  Center  (SRDC)  at  the  Naval  Postgraduate 
School  (NPS)  has  a  high  quality  adaptive  optics  testbed,  which  was  developed  in 
conjunction  with  personnel  from  the  Naval  Research  Laboratory  (NRL)  in  Albuquerque, 
NM.  The  purpose  of  the  testbed  is  to  apply  control  techniques  to  remove  aberrations 
from  an  optical  wavefront.  One  application  that  demonstrates  this  purpose  is  an  imagery 
application,  where  adaptive  optics  are  used  to  improve  the  image  quality  of  an  object  of 
interest.  Figure  6  and  Figure  7  show  an  example  of  the  improvement  in  image  quality 
that  can  be  obtained  using  this  technique.  The  four  lights  at  the  bottom  are  LEDs  used  to 
illuminate  the  object. 
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Figure  7  Corrected  AO  Testbed  Image. 


The  primary  components  of  the  testbed  include  two  deformable  mirrors  and  two 

Shack-FIartmann  wavefront  sensors,  two  fast  steering  mirrors,  a  reference  laser,  a  source 

object,  and  a  science  camera.  The  reference  and  object  light  beams  are  conveyed  around 
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the  table  by  several  lenses,  mirrors,  and  beam  splitters  which  relay  the  pupil  plane  to  the 
different  components.  A  schematic  of  the  testbed  is  shown  in  Figure  8  and  a  picture  is 
shown  in  Figure  9. 


Light  from  the  reference  laser  beam  and  the  object  light  are  combined  at  the  first 
beam  splitter.  This  combined  light  is  then  reflected  off  a  Micromachined  Deformable 
Mirror  (MMDM).  A  Shack-Hartmann  wavefront  sensor  measures  a  portion  of  the 
reflected  reference  light  source  and  the  rest  continues  through  the  system.  The  Shack- 
Hartmann  wavefront  sensor  sends  its  CMOS  camera  image  to  the  control  computer, 
which  converts  the  CMOS  image  into  a  set  of  wavefront  slopes.  A  control  algorithm  is 
used  to  determine  the  correction  necessary  to  make  the  wavefront  planar,  and  a  set  of 
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voltages  is  sent  to  the  deformable  mirror  to  implement  the  correction.  This  is  the  first 
closed  loop  control  system  on  the  testbed. 

A  filter  in  front  of  the  Shack-Hartmann  wavefront  sensor  removes  the  white 
object  light  prior  to  hitting  the  lenslet  array.  Therefore,  any  wavefront  aberrations  are 
measured  on  the  reference  beam,  not  the  object  light.  There  are  two  important 
assumptions  of  this  setup.  The  first  is  that  the  same  aberrations  are  present  in  the 
reference  beam  as  the  image  beam.  Second,  correcting  the  reference  beam  will  also 
correct  the  object  light.  It  is  possible  that  since  the  laser  light  is  monochromatic  while  the 
object  light  contains  all  visible  wavelengths,  some  of  the  optical  elements  on  the  table 
could  cause  chromatic  aberrations  in  the  object  light  that  are  not  present  in  the  laser 
beam.  These  differences  are  assumed  to  be  negligible. 

After  the  first  closed  loop  wavefront  correction  process,  the  combined  light  beam 
next  encounters  two  Fast  Steering  Mirrors  (FSMs).  Including  just  one  FSM  allows  tip/tilt 
jitter  aberrations  to  be  removed.  The  presence  of  two  FSMs  allows  the  injection  of  a 
disturbance  with  the  first  FSM  to  simulate  an  on-board  jitter,  and  the  removal  of  the  jitter 
with  the  second  FSM.  Of  course,  in  their  nominal  configuration,  when  no  commands  are 
sent  to  the  actuators  on  the  FSMs,  they  behave  as  ordinary  planar  mirrors.  The  Position 
Sensing  Detector  (PSD)  in  the  schematic  is  the  wavefront  sensor  used  to  close  the  control 
loop  with  the  second  FSM. 

The  combined  light  beam  next  encounters  a  second  deformable  mirror,  this  one 
having  piezoelectric  actuators  (PDM).  Like  the  first  DM,  the  light  reflected  off  it  is 
imaged  by  a  Shack-Hartmann  wavefront  sensor  and  analyzed  by  a  second  control 
computer  to  detennine  the  proper  correction  to  be  applied  by  the  DM.  This  completes  the 
final  closed  loop  system  on  the  testbed.  Finally,  the  reference  laser  beam  is  filtered  out 
and  a  science  camera  is  able  to  take  an  image  of  the  object  light  as  seen  in  Figure  7. 

2.  Hardware 

This  section  will  give  specifics  on  some  of  the  AO  testbed  hardware.  The  control 
laws  developed  involve  only  the  components  in  the  first  closed  loop  system;  the  MMDM 
and  one  Shack-Hartmann  wavefront  sensor.  The  other  components  did  not  have  an  active 
role  in  verifying  the  performance  of  the  control  laws. 
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a.  Laser 

The  laser  on  the  testbed  is  a  red  HeNe  laser  with  a  wavelength  of  632.8 
mn.  It  is  manufactured  by  JDS  Uniphase  and  the  model  number  is  1 137P. 

b.  Micromachined  Deformable  Mirror 

The  MMDM  is  manufactured  by  OKO  Technologies  of  Delft,  the 
Netherlands.  The  following  information  is  taken  from  the  user  manual  [19].  The 
aperture  is  15  mm  in  diameter.  Thirty-seven  actuators  are  located  under  a  12  mm 
diameter  of  the  membrane  with  1.8  mm  center-to-center  spacing  between  actuators.  A 
diagram  is  shown  in  Figure  10.  The  mirror  itself  is  a  silicon  chip  coated  with  a  silicon 
nitride  membrane.  The  numbering  and  two-dimensional  positioning  of  the  actuators  is 
shown  in  Figure  1 1 . 


Figure  10  MMDM  Schematic  Section  From  [19]. 
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Figure  1 1  MMDM  Actuator  Locations  and  Positioning  From  [19]. 


c.  Shack-Hartmann  Wavefront  Sensor 

OKO  Technologies  [19]  also  provides  the  Shack-Hartmann  wavefront 
sensor  on  the  testbed.  It  is  a  14- inch  CMOS  A601f  camera  manufactured  by  the  German 
company  Basler  with  a  Hartmann  lenslet  mask  from  OKO  installed.  The  sensor  has  656 
x  491  pixels  and  a  frame  rate  of  30  frames  per  second.  The  lenslet  mask  has  127  lenslets 
laid  out  in  a  hexagonal  array  and  coated  with  a  fused  silica  coating.  The  aperture  size  is 
3.5  mm  and  the  subaperture  of  each  lenslet  is  100  pm.  The  appearance  of  the  array  is 
similar  to  that  of  one  of  the  segments  in  the  segmented  space  telescope  model  array 
shown  in  Figure  5,  only  with  127  lenslets  vs.  61. 

d.  Control  Computer 

The  control  algorithms  are  implemented  on  a  standard  desktop  PC  and 
Matlab  software.  The  computer  is  a  Pentium  D  with  a  2.8  GHz  processor  and  1.0  GB  of 
RAM.  The  operating  system  is  Windows  XP  version  2002,  service  pack  2.  The  Matlab 
version  is  R2007a.  The  Basler  camera  is  connected  via  fire  wire  to  the  computer,  and 
two  PCI  cards  connect  the  computer  to  the  MMDM.  One  card  controls  19  actuators,  the 
other,  18. 
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III.  ROBUST  CONTROL 


A.  TRADITIONAL  CONTROL  APPROACH 

The  purpose  of  any  control  system  is  to  provide  stability  and  achieve  certain 
performance  objectives.  For  adaptive  optics  systems,  classical  control  techniques  have 
always  been  adequate.  Their  simpler  design  process  and  more  intuitive  nature  make  them 
the  controller  of  choice  for  most  applications.  The  following  will  give  a  brief  overview 
of  the  traditional  AO  system  control  approach  and  identify  the  shortcomings  that  led  to 
this  research  into  more  advanced  control  techniques. 

Traditionally,  the  adaptive  optic  in  large  ground  based  AO  systems  is  treated  as  a 
static  system  with  no  dynamics.  In  the  case  of  a  deformable  mirror,  it  is  assumed  that  the 
coupling  between  actuators  is  static.  Classical  control  approaches  neglect  the  dynamics 
associated  with  this  coupling  thereby  simplifying  the  control  process  [3],  [5].  Obviously, 
as  the  optics  for  space  telescopes  become  lighter  as  discussed  in  Chapter  I,  and  as  more 
actuators  are  packed  into  a  given  area  to  provide  higher  spatial  correction,  this 
simplifying  assumption  becomes  less  realistic.  For  systems  with  a  high  degree  of  internal 
dynamics,  instability  can  result  from  a  controller  designed  from  a  static  model. 

At  the  basis  of  the  control  algorithm,  the  reconstruction,  or  influence  matrix, 
relates  the  actuators  with  the  aberrations  in  the  wavefront.  The  reconstruction  matrix  is 
formed  by  starting  with  the  mirror  in  a  neutral  position  and  then  actuating  or  poking  each 
actuator  in  turn  while  measuring  the  effect  on  the  wavefront.  For  this  reason,  the 
reconstruction  matrix  is  often  referred  to  as  the  poke  matrix.  As  an  example,  consider  a 
deformable  mirror  with  10  actuators  and  a  Shack-Hartmann  wavefront  sensor  with  20 
lenslets.  Assume  the  neutral  position  of  each  actuator  is  represented  by  a  0,  and  it  can  be 
commanded  to  any  value  from  -1  to  +1.  Actuator  1  is  given  full  positive  stroke  and  an 
image  is  taken  with  the  Shack-Hartmann  wavefront  sensor.  This  measurement  yields  a 
vector  of  40  measurements  (20  x  and  20  y.)  Actuator  1  is  reset  to  0,  and  the  process  is 
repeated  for  actuator  2.  The  end  result  is  a  reconstruction  matrix  with  dimensions  of  [40 
x  10]  where  each  column  represents  the  x  and  y  pixel  displacements  on  the  camera  sensor 
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for  each  individual  actuator.  These  x  and  y  pixel  displacements  can  be  easily  converted 
to  x  and  y  slopes  by  multiplying  by  the  physical  pixel  width  and  dividing  by  the  camera 
focal  length. 


This  process  creates  a  map  of  how  an  actuator  position  or  command  affects  the 
wavefront,  but  it  neglects  the  structural  dynamics  of  the  mirror.  The  mathematical 
expression  of  the  reconstruction  matrix,  R ,  is  given  as 
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where  S  is  the  slope  in  the  x  or  y  direction,  n  is  the  number  of  lenslets, 
number  of  actuators.  The  relationship  between  the  actuator  input  vector,  u , 
measurement  vector,  S ,  is  given  by 


(3.1) 


and  m  is  the 
and  the  slope 


S  =  Ru  (3.2) 

Since  the  slope  is  the  measured  quantity,  the  pseudo-inverse  of  the  reconstruction  matrix 
is  used  to  determine  the  control  inputs.  One  possible  control  approach  is  to  subtract  the 
product  of  the  inverted  reconstruction  matrix  (found  with  a  pseudo-inverse  technique) 
and  measured  slope  from  the  previous  control  input  as  shown  in  Equation  (3.3).  This  will 
drive  the  control  signal  to  a  steady  state  value  and  the  wavefront  slope  to  zero,  resulting 
in  a  planar  wavefront. 

a m=^u-R's  (3.3) 

Variations  on  how  the  wavefront  is  detennined  from  the  slopes  (e.g.,  modal  or 
zonal)  and  variations  on  the  feedback  method  (e.g.,  direct,  indirect,  iterative,  etc.)  can  be 
employed  to  improve  the  performance  for  this  traditional  control  law  [3], 


As  stated  before,  the  classical  control  approach  ignores  the  structural  dynamics  of 
a  system.  Even  if  a  more  complex  multi-input,  multi-output  control  approach  was  used 
that  accounted  for  dynamics  and  coupling  between  actuators,  it  would  still  need  to 
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perform  satisfactorily  in  the  presence  of  external  disturbances  and  unmodeled  dynamics. 
These  uncertainties  give  credence  to  the  idea  of  designing  robust  controllers  for  AO 
systems. 

B.  JUSTIFICATION  FOR  ROBUST  CONTROL 

Robust  control  addresses  the  problem  of  designing  a  control  system  that  yields  a 
desired  perfonnance  in  the  presence  of  model  uncertainties  and  external  disturbances. 
For  the  purposes  of  this  research,  the  plant  to  be  controlled  is  a  state  space  model  of  the 
following  form 

\x  =  Ax  +  Bu 

G  =  \  n  n  (3.4) 

[  y  =  Cx  +  Du 

where  x  is  the  state,  y  is  the  measured  output,  and  a  is  the  control  input.  The  robust 
control  synthesis  process  augments  the  plant,  G,  into  a  modified  system,  P,  with  two 
input  vectors  and  two  output  vectors.  The  process  of  constructing  the  modified  plant,  P, 
will  be  described  later.  The  robust  controller,  K,  is  then  designed  for  this  augmented 
plant  shown  in  Figure  12 


Figure  12  Standard  Robust  Control  System. 


where  the  plant,  P,  and  controller,  K,  are  both  state  space  models,  r  is  the  external  input 
(which  contains  a  reference  input,  external  disturbances,  and  noise),  z  is  the  performance 
measure,  u  is  the  control  input,  and  e  is  the  error  between  the  reference  input  and 
measured  output.  For  the  specific  case  of  the  segmented  space  telescope  model  used  for 
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this  research,  the  state  space  plant  model  describes  both  the  actuators  and  the  Shack- 
Hartmann  wavefront  sensor  outputs.  The  external  inputs,  r,  are  all  zero.  The  reference 
input  is  usually  zero  for  adaptive  optics  applications  since  a  planar  wavefront  is  desired; 
however,  in  some  situations  a  non-zero  reference  input  may  be  used  if  a  particular 
aberration  mode  is  desired.  Since  the  error  signal,  e,  is  the  difference  between  the 
measured  output,  y,  and  the  reference,  r,  and  the  reference  is  zero,  the  error  becomes 
simply  the  negative  of  the  measured  output.  The  control  input,  u,  is  a  vector  of  actuator 
voltages.  For  this  research  the  performance  measure,  z,  is  comprised  of  two  components, 
the  first  is  the  error  signal,  e,  and  the  second  is  the  measured  output,  y.  These  two 
performance  measures  will  be  penalized  or  weighted  differently  for  the  controller 
synthesis  process.  This  weighting  process  will  be  discussed  shortly.  The  end  goal  then 
of  the  robust  controller  synthesis  process  is  to  minimize  the  oo-nonn  of  the  operators 
from  r  — >  z  . 

The  overall  goal  of  a  robust  controller  is  to  provide  stability  for  a  system  in  the 
presence  of  uncertainty  or  unmodeled  dynamics  as  discussed  in  Chapter  I.  Ideally,  a 
robust  controller  will  be  able  to  handle  these  uncertainties  better  than  other  controllers 
such  as  classical  PID  controllers  or  optimal  linear  quadratic  regulators.  When  speaking 
of  robustness,  there  are  two  different  perfonnance  criteria  implied:  the  ability  to  handle 
uncertainty  in  the  model,  and  the  ability  to  reject  disturbances  and  noise.  What  separates 
the  robust  control  design  process  from  a  traditional  control  design  process  is  that  these 
uncertainties  are  explicitly  accounted  for  in  the  design  process  to  improve  performance  in 
the  face  of  uncertainties  that  are  poorly  understood  or  modeled.  The  //y  robust  control 
method  minimizes  the  upper  bound  of  the  perfonnance  signal. 

What  follows  are  some  of  the  mathematical  foundations  that  form  the  backbone  of 
robust  control  theory.  Much  of  the  material  is  general  in  nature  and  not  specific  to  the 
research  performed  on  the  analytical  segmented  space  telescope  model.  The  topics 
include  system  norms,  sensitivity  and  complimentary  sensitivity  functions,  uncertainty 
modeling,  and  some  of  the  tools  necessary  for  the  fonnulation  of  the  appropriate 
augmented  system  matrix,  P,  to  which  the  controller  is  synthesized.  The  application  of 

these  techniques  to  the  analytical  model  will  be  covered  in  Chapter  V. 
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C.  MATHEMATICAL  FOUNDATIONS 
1.  Norms 


The  basic  premise  of  robust  control  involves  minimizing  a  system  nonn.  The 
type  of  norm  is  dependent  on  the  type  of  robust  controller  being  designed.  A  general 
overview  of  norms  follows.  From  a  conceptual  point  of  view,  a  norm  is  nothing  more 
than  a  measure  of  how  large  something  is,  such  as  the  size  of  a  vector.  The  general  form 
of  a  vector  norm  in  □  p  as  shown  in  [16]  is  given  as 


f  V 

II 


\Vp 


1  <  p  <  oo 


(3.5) 


V  i= 1  J 

where  p  takes  on  the  value  of  the  desired  norm.  The  1-norm  and  2-norm  (also  known  as 
the  Euclidean  norm)  are  common  norms.  In  the  case  of  robust  control,  a  common  nonn, 
and  the  one  used  in  this  research,  is  the  infinity-norm,  given  as 

llxll  :=  max  lx,.  I  (3.6) 

II  llco  \<i<tt  !| 


The  standard  designation  for  this  type  of  controller  is  //y ,  where  the  H  refers  to  the 


Hardy  space  for  stable  systems. 


For  matrix  operations,  the  induced  nonn  becomes  more  useful,  where  the  norm 
can  alternately  be  thought  of  as  a  gain.  Given  the  following  simple  transformation  where 
A  is  the  gain  or  amplification  of  the  system 

y  =  Ax  (3.7) 

the  induced  matrix  nonn  is  given  as 

1 1 A  JC 1 1 

||.<4||  =  max  p  (3.8) 

p  x*°  x 


2.  Singular  Values 

Singular  values  are  used  extensively  for  robust  control  synthesis.  Using  a 
standard  definition  for  linear  algebra,  singular  values  are  the  square  root  of  the 
eigenvalues  of  a  matrix  times  its  adjoint  (i.e.,  complex  conjugate  transpose)  as  shown  by 

ai  =  ^Ai(A*A)  (3.9) 
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It  can  be  shown  that  the  largest  gain  for  any  input  direction  is  equal  to  the 
maximum  singular  value,  cr,  and  that  the  smallest  gain  for  any  input  is  equal  to  the 
minimum  singular  value,  a .  The  Hx  norm  can  also  be  defined  as  the  least  upper  bound 
on  the  maximum  singular  value.  The  Hx  norm  for  a  linear  time  invariant  system,  G,  is 
given  as 

||G|LDsupa(G(»)  (3.10) 

CO 

For  multi-input,  multi-output  systems,  a  plot  of  the  maximum  singular  values  is 
analogous  to  the  Bode  magnitude  plot  for  single-input,  single-output  systems. 

3.  Small  Gain  Theorem 

One  of  the  most  fundamental  concepts  for  robust  control  is  the  small  gain 
theorem,  which  follows  from  the  above  discussion  of  norms.  For  stability  analysis,  the 
small  gain  theorem  replaces  the  traditional  concepts  of  gain  and  phase  margin.  Consider 
the  system  in  Figure  12  with  uncertain  elements  A  added  as  shown  in  Figure  13. 


T 


Figure  13  System  with  Uncertainty. 
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Let  the  region  inside  the  dashed  line,  which  is  identical  to  the  system  in  Figure  12, 
be  represented  as  M.  The  small  gain  theorem  requires  than  the  gain  (norm)  of  the 
uncertainty,  A,  multiplied  by  the  gain  of  the  system,  M,  must  be  less  than  1  for  the 
system  to  be  robustly  stable  as  shown  by 

||a||*||m||<i  va  (3.11) 

where  the  *  represents  a  multiplication. 

Another  common  form  of  the  small  gain  theorem  is  given  as 

HLS  yr  W  \Ml<r  (3.12) 

where  y  is  a  bound  and  y  >  0 . 

4.  Sensitivity  Functions 

Two  other  necessary  concepts  for  robust  control  design  are  the  sensitivity 
function  and  the  complementary  sensitivity  function.  These  functions  allow  the  designer 
to  shape  the  frequency  response  of  the  system  in  a  manner  similar  to  classical 
loopshaping  techniques.  The  sensitivity  function  is  important  for  shaping  the  response  of 
lower  frequencies,  while  the  complimentary  sensitivity  function  helps  shape  the  high 
frequency  response  [20].  Mathematically,  for  a  system  consisting  of  a  plant,  G,  and  a 
controller  ,K,  with  negative  feedback,  the  sensitivity  function,  S,  has  the  following  form 

S  =  (I  +  GK r1  (3.13) 

From  this,  it  follows  that  the  complementary  sensitivity  function,  T,  has  the  following 
form 

T  =  GK(I  +  GK)1  (3.14) 

Due  to  their  complimentary  nature,  the  following  is  true 

S  +  T  =  1  (3.15) 

A  control  design  process  that  incorporates  loopshaping  with  both  the  sensitivity 
and  complementary  sensitivity  functions  is  often  referred  to  as  a  mixed-sensitivity 
problem  [20].  In  this  instance,  the  sensitivity  function  influences  the  performance  of  the 
controller,  and  the  complementary  sensitivity  function  influences  the  stability  of  the 
system  at  higher  frequencies.  Figure  14  gives  a  representation  of  how  loopshaping  of  the 
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sensitivity  and  complementary  sensitivity  functions  can  affect  the  controller  /  plant 
system,  represented  by  L  ( L  =  GK )  in  the  figure. 


Performance 


Robustness 


Figure  14  Robust  Control  Loopshaping  From  [20]. 


At  low  frequencies,  the  loop  transfer  function,  L,  is  weighted  to  follow  the  inverse 
of  the  sensitivity  function.  The  sensitivity  function  of  the  system  describes  how  well  the 
controller  does  at  providing  for  disturbance  rejection,  tracking  errors,  and  negating  any 
plant  parametric  variations.  Shaping  the  sensitivity  function  indirectly  shapes  the  loop 
transfer  function,  and  in  some  cases  vice  versa.  At  higher  frequencies,  the  shape  of  the 
complementary  sensitivity  function  describes  the  transient  response  and  stability.  The 
complementary  sensitivity  function  is  shaped  by  weighting  functions  in  the  same  way  as 
the  sensitivity  function. 

5.  Weighting  Functions 

At  this  point,  it  is  necessary  to  introduce  a  set  of  weighting  functions  that  will 
manipulate  or  shape  the  sensitivity  and  complementary  sensitivity  functions.  The 
weighting  functions  can  take  whatever  fonn  is  appropriate  for  the  system,  but  for  the 
purpose  of  this  research,  first  order  transfer  functions  are  satisfactory.  The  first 
weighting  function,  introduced  as  Wj,  is  a  penalty  on  the  error  and  affects  the  lower 
frequencies  and  shape  of  the  sensitivity  function.  A  second  weighting  function,  W2,  is  a 
penalty  on  the  control  signal  but  is  not  used  in  this  research.  The  final  weighting 
function,  Wj,  is  a  penalty  on  the  measured  output  and  affects  the  complementary 
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sensitivity  function  and  higher  frequencies.  Figure  15  shows  a  schematic  of  a  nominal 
system  with  a  plant,  a  controller,  and  weights  Wj  and  W3.  Note  that  the  transfer  functions 
shown  in  the  figure  merely  represents  a  generic  transfer  function  and  not  the  actual 
weights  used. 


Figure  15  Schematic  of  System  with  Weighting  Functions. 

The  weighted  error  and  measured  output  signals  together  form  the  performance 
measure  previously  identified  in  Figure  12  as  z.  The  error  that  is  weighted  by  Wx  is 
minimized  in  order  to  achieve  the  desired  performance  by  influencing  the  sensitivity 
function.  The  output  that  is  weighted  by  W3  is  minimized  in  order  to  achieve  stability  in 

the  presence  of  higher  frequency  disturbances  by  influencing  the  complementary 
sensitivity  function.  Satisfying  both  of  these  goals  simultaneously  may  not  be  possible; 
thus  an  iterative  approach  to  synthesizing  the  robust  controller  is  used.  If  the  weights 
selected  are  too  restrictive  on  the  design,  a  robust  controller  meeting  the  design 
requirements  cannot  be  synthesized. 

The  closed  loop  system  with  weight  Wj  is  said  to  have  robust  performance  if  the 
following  condition  of  the  small  gain  theorem  is  met 

|^|<|i  +  ga:|  v®  (3.16) 

Similarly,  when  considering  W3  and  the  complementary  sensitivity  function,  the  system  is 
robustly  stable  if  the  following  condition  is  true 

\W3GK\<\l  +  GK\  Vfl)  (3.17) 
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Combing  equations  (3.16)  and  (3.17),  the  system  has  both  robust  performance  and 
stability  if  the  following  is  true 

l^l+l^rllL  <1  (3.18) 

Designs  that  employ  a  weight  on  both  the  sensitivity  function  and  the 
complementary  sensitivity  function  are  referred  to  as  weighted  mixed  sensitivity 
problems  [20]. 

D.  MODEL  UNCERTAINTY 

The  main  advantage  of  designing  a  robust  controller  over  other  simpler  classical 
techniques  is  the  ability  to  handle  model  uncertainties.  While  certain  designs  allow  for 
wide  latitude  in  handling  uncertainties,  it  would  be  nearly  impossible  to  design  one 
controller  for  all  possible  uncertainties  a  system  might  encounter  over  its  lifetime. 
Therefore,  it  is  necessary  to  understand  what  kinds  of  uncertainty  are  likely,  where  they 
come  from,  and  how  to  bound  the  uncertainties  to  simplify  some  of  the  design 
constraints.  Uncertainty  can  be  added  to  a  system  in  two  different  ways,  an  additive 
fashion  or  a  multiplicative  fashion.  The  following  two  figures  depict  these  two  methods. 
Perturbations  of  these  types  can  be  added  to  any  combination  of  inputs  or  outputs  to  the 
system. 


ga  =  Go  +  aa 


Figure  16  Additive  Uncertainty  From  [20]. 


Am 

G„ 

■  + 

+  V- 

ga  -  Go(I  +  Am) 


Figure  17  Multiplicative  Uncertainty  From  [20]. 
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The  first  consideration  in  defining  the  uncertainty  to  be  added  to  a  model  is  to 
define  it  as  structured  or  unstructured.  An  unstructured  uncertainty  is  a  bounded 
perturbation  that  is  included  or  added  to  the  model,  but  where  very  little  of  the  nature  of 
the  perturbation  is  known.  For  example,  there  is  no  transfer  function  that  describes  the 
general  behavior  of  an  unstructured  uncertainty. 

Structured  uncertainties  involve  perturbations  that  are  more  understood  than 
unstructured.  In  general,  there  are  more  constraints  on  the  nature  of  a  structured 
uncertainty  than  an  unstructured  uncertainty.  Structured  uncertainties  can  result  from 
bounded  variations  in  model  parameters  such  as  changes  in  the  physical  components  of  a 
system  over  time  due  to  wear  and  tear.  They  also  result  from  combinations  of  multiple 
uncertainties,  or  both  additive  and  multiplicative  perturbations.  Figure  18  shows  a 
diagram  of  typical  structured  and  unstructured  model  uncertainties  [17]. 


•  U:  Unstructured 

•  S:  Structured 

Figure  18  Structured  and  Unstructured  Model  Uncertainty  Sources  From  [20]. 
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In  addition  to  helping  shape  the  sensitivity  and  complementary  sensitivity 
functions,  the  weighting  functions  also  bound  the  system  uncertainty  [20],  For  this 
reason,  proper  selection  of  the  weighting  functions  becomes  a  crucial  part  of  robust 
control  design. 


The  segmented  space  telescope  model  used  in  this  research  does  not  contain  any 
unstructured  uncertainties.  On  the  other  hand,  uncertainties  on  some  of  the  physical 
parameters  such  as  modal  frequencies  and  damping  are  assumed.  This  affects  the  lower 
left  quadrant  where  the  squares  of  the  natural  frequencies  are  found  in  the  following 
system  A  matrix 


A  = 


X 

0 

/ 

"2  fan 

'■) 

(3.19) 


Uncertainty  in  the  natural  frequencies  is  a  logical  concern  for  a  lightweight  space 
telescope  model  such  as  this.  Therefore,  a  percentage  of  uncertainty  can  be  added  to  the 
natural  frequencies  to  allow  for  explicit  design  of  a  robust  controller  based  on  an 
uncertain  model. 


E.  CONTROLLER  SYNTHESIS 

The  robust  controllers  designed  for  this  research  were  based  on  minimizing  the 
closed  loop  system  oo  -norm.  Other  robust  control  synthesis  methods  exist  such  as  those 
based  on  a  2-norm,  but  they  will  not  be  discussed  here.  The  go  -nonn  method  used  here  is 
a  direct  method  denoted  as  Hx .  This  method  designs  a  controller  based  on  a  nominal 
augmented  plant  as  previously  shown  in  Figure  12. 

The  Hx  robust  control  design  technique  requires  a  properly  constructed  plant 

model  to  synthesize  a  robust  controller.  In  order  to  achieve  the  robust  performance  and 
stability  design  objectives,  the  plant  must  be  augmented  with  the  appropriate  weighting 
functions.  A  representation  of  this  augmented  nominal  plant  is  shown  in  Figure  19. 
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Input  2  Output  2 


Figure  19  Augmented  Nominal  2-Input,  2-Output  Plant. 


This  augmented  plant  is  often  referred  to  as  the  2-input,  2-output  model.  The  two 
vector  inputs  to  the  augmented  plant  are  the  reference  input  and  the  control  input,  which 
the  controller  to  be  designed  will  provide.  One  of  the  vector  outputs  is  the  weighted  error 
and  weighted  measured  output  signal  combined,  which  is  otherwise  known  as  the 
performance  measure.  The  other  vector  output  is  the  error  signal  that  is  used  as  the  input 
for  the  controller. 


Once  this  2-input,  2-output  model  has  been  constructed,  the  goal  is  to  minimize 
the  oo  -nonn  of  the  error  signal  to  the  command  or  disturbance  input.  A  representation  of 
how  the  weighting  functions  are  included  in  the  augmented  model,  P,  is  given  as 


P  = 


W1 

-WIG 

0 

W3G 

I 

-G 

(3.20) 


Since  the  original  plant,  G,  was  a  state  space  model,  P  is  also  a  state  space  model 
and  can  be  partitioned  as 


(A 

B i 

b2 

P  = 

Cy 

Ai 

Dn 

[c2 

D2\ 

A2 

(3.21) 


where  the  A,  B,  C,  and  D  matrices  are  the  augmented  state  space  matrices.  The 
controllability  and  observability  of  the  system  can  be  determined  using  the  submatrices 
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Bx,  B2,  etc.  As  shown  in  [16],  the  following  assumptions  on  the  partitions  of  P  are  made 
before  proceeding  with  the  controller  synthesis 

•  (A,B j )  is  controllable 

•  (C15H)  is  observable 

•  (A,B2)  is  stabilizable 

•  (C2,A)  is  detectable 

•  D’AC,  Dn)  =  (  0  I) 


(Bx  ^ 

fo" 

vA  2  j 

Dh  = 

Ia 

An  admissible  controller  exists  if  and  only  if  the  following  conditions  are  met: 

•  He  ( Ric )  and  X  :=  Ric(H)  >  0 

•  J  e  (Ric)  and  Y  :=  Ric(J)  >  0 

•  p(XY)<Y 2 

where  H  and  J  are  the  following  Hamiltonian  matrices 

(  A  r2BX-B2BT2] 

{ <CX  -AT  J 


(3.22) 


J  := 


At 

k-bX 


T ^  \ 
2 


r^qq-c'c 

-A 


(3.23) 


y  is  a  performance  bound  with  y  >  0 ,  Ric  is  the  solution  to  an  algebraic  Riccati  equation 
of  the  form 


ATX  +  XA  +  XBX  +  C  =  0  (3.24) 

and  p  is  the  spectral  radius  with 

p(A):=  max  U.  I  (3.25) 

v  7  1  <i<n  '  1 

The  robust  control  design  algorithm  iterates  on  the  value  of  y  between  0  and 
infinity  with  a  bisection  algorithm  until  such  time  as  a  solution  to  the  Riccati  equation 
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cannot  be  found.  Once  the  specified  tolerance  threshold  has  been  reached  the  iteration 
stops  and  the  controller,  K,  can  be  computed  from  the  following  equations 


where 


f  Af 

-ZL 

zb2 

K  = 

M 

0 

I 

-C 

\  ^2 

I 

0 

M  =  -B‘  X 


L  =  -YCl 


Z  =  [l-y~2YX)1 


(3.26) 


(3.27) 

(3.28) 

(3.29) 


Af  =A  +  y-2BxB{X  +  B2M  +  ZLC2  (3.30) 

Figure  20  shows  the  synthesized  controller,  K,  attached  to  the  augmented  plant 
from  Figure  19.  Figure  20  is  now  identical  to  the  original  system  presented  in  Figure  12. 


Figure  20  Augmented  Plant  and  Controller. 

The  controller  synthesis  consists  of  designing  the  controller,  K,  such  that  the 
following  cost  function  is  minimized 
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\Kl  <r  (3.3i) 

where  T  is  the  closed  loop  transfer  function.  After  this  design  process  has  generated  the 
Hx  controller,  it  can  be  implemented  as  a  continuous  system,  or  discretized  prior  to 

implementation.  The  application  of  these  techniques  to  the  segmented  space  telescope 
model  will  be  discussed  in  Chapter  V. 


42 


IV.  MODEL  REDUCTION  TECHNIQUES 


A.  INPUT/OUTPUT  REDUCTION 

The  robust  controller  synthesis  process  described  in  Chapter  III  can  be  rather 
computationally  prohibitive.  In  the  case  of  the  analytical  segmented  space  telescope 
model,  the  332  states,  997  inputs,  and  936  outputs  are  more  than  an  average  desktop 
computer  can  process.  Therefore,  some  kind  of  model  reduction  is  necessary.  This 
chapter  will  present  an  innovative  method  developed  to  reduce  the  inputs  and  outputs  as 
well  as  a  standard  method  of  state  reduction. 

1.  Singular  Value  Decomposition 

The  principle  behind  the  reduction  method  developed  is  the  projection  of  a  vector 
into  a  subspace.  In  this  case,  a  large  vector  of  input  values  is  projected  into  a  vector  of 
many  fewer  values.  If  the  projection  is  chosen  carefully,  then  the  important  or  significant 
values  of  the  original  vector  are  preserved  so  that  when  the  reduced  vector  is  expanded 
back  to  the  original  number  of  values,  it  is  still  representative.  The  Singular  Value 
Decomposition  (SVD)  is  one  such  technique  that,  by  its  very  nature,  preserves  the 
significant  values. 

In  general,  the  SVD  has  the  following  form 

SVD(H)  =  UT.Vr  (4.1) 

where  H,  U,  E ,  and  V  are  all  matrices,  and  I  is  a  diagonal  matrix  with  the  singular 
values  arranged  in  decreasing  size.  The  matrices  U  and  V  are  also  ordered  such  that  the 
vector  associated  with  the  largest  singular  value  is  in  column  one  and  so  on.  A  plot  of  all 
the  singular  values  helps  to  identify  how  many  singular  values  have  a  sufficient  impact 
on  the  final  system  so  that  a  logical  place  to  truncate  the  singular  value  matrix  can  be 
found.  For  example,  if  there  are  40  total  singular  values,  but  the  last  25  are  significantly 
smaller  than  the  first  15,  those  last  25  values  can  be  truncated  and  the  most  important 
aspects  of  the  original  system  are  still  preserved.  Since  the  column  vectors  of  the  U  and 
V  matrices  are  also  rank  ordered  into  decreasing  importance,  only  the  first  15  column 
vectors  of  those  matrices  will  be  retained.  As  will  be  shown  in  Chapter  V,  when  H  is 
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defined  to  be  the  transfer  function  of  the  inputs  to  the  outputs  for  the  analytical  space 
telescope  model,  the  large  number  of  inputs  or  outputs  of  the  original  system  can  be 
projected  into  the  column  vectors  of  the  truncated  U  and  V  matrices.  The  inputs  are 
projected  via  the  V matrix  and  the  outputs  are  projected  via  the  U  matrix. 

Once  a  controller  has  been  synthesized  for  this  reduced  input/output  model,  the 
controller  outputs  can  be  projected  back  into  the  total  number  of  original  plant  inputs 
through  the  truncated  V matrix.  Similarly,  the  original  plant  outputs  are  projected  into  the 
truncated  U  matrix  so  that  they  may  be  inputs  for  the  reduced  controller. 


2.  Zernike  Polynomials 


An  innovative  approach  developed  in  this  dissertation  is  the  use  of  Zemike 
polynomials  for  model  output  reduction.  Zemike  polynomials  are  orthogonal  and  nonnal 
about  a  unit  circle,  making  them  especially  applicable  to  optical  applications  where  the 
optical  components  have  circular  apertures.  Their  orthogonal  and  normal  nature  as 
described  in  [21]  is  shown  by 

j  K(p)K’(p)pdp  =  7<S„,  (4.2) 


where  R"'  is  the  radial  polynomial  and  8  is  the  Kronecker  delta.  It  is  worth  noting, 


however,  that  Zemike  polynomials  are  but  one  of  an  infinite  number  of  polynomial  sets 
that  are  nonnal  about  a  unit  circle.  Furthermore,  variations  on  the  Zernike  polynomials 
exist  depending  on  their  particular  application  and  usage.  Derivations  of  the  standard 
form  of  Zemike  polynomials  can  be  found  in  [22],  [23],  A  slightly  modified  form  of 
Zernike  polynomials  known  as  the  fringe  Zernikes  was  developed  at  the  University  of 
Arizona  and  is  used  for  this  research  [24].  The  main  difference  is  that  the  fringe  Zernikes 
have  been  normalized  to  unity  magnitude  at  the  edge  of  the  pupil.  The  expansion  series 
of  the  fringe  Zernike  polynomials  is  given  by 

z:(p,e)=K(py“  (4.3) 

where 


’lj^  (2  n-m—s)l  ,,  , 

K(p)  =  Z(-i)  {  w — 


5=0 


s\{n -x)!(n  -m  -x)! ' 


(4.4) 
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where  Z”!  and  R"‘  both  represent  the  Zernike  polynomial,  n  and  m  are  both  integers 
with/?  >  m  ,  ?7  >  0  ,  and  p  is  a  function  of  the  aperture  radii.  The  index,  n,  defines  the 
order  of  the  radial  power,  and  m  is  the  meridional  frequency.  The  meridional  frequency 
is  in  reference  to  the  meridional  or  tangential  plane  of  the  aperture,  which  is  a  plane  that 
contains  both  the  optical  axis  and  the  chief  ray.  The  chief  ray  is  the  ray  that  passes 
through  the  center  of  the  aperture.  Due  to  the  alternating  cosine  and  sine  terms  in  the 
expansion,  the  polynomials  can  be  defined  somewhat  differently  based  on  how  the  index 
m  is  handled.  Figure  21  shows  a  pictorial  representation  of  the  first  20  Zemike  modes 
using  the  standard  fonn. 


Figure  2 1  Twenty  Common  Zernike  Modes  From  [25]. 


Zernike  polynomials  are  usually  derived  in  polar  form  but  can  be  converted  to  the 
Cartesian  form  used  in  this  research  through  the  following  equations: 


45 


x  =  /?cos(#) 
v  =  /?sin(<9) 


(4.5) 


(4.6) 


Both  the  standard  fonn  and  the  fringe  Zemikes  have  the  same  orthogonality 
properties.  The  main  factor  that  sets  the  Zernike  polynomials  apart  from  other 
polynomial  sets  is  a  group  of  invariance  properties  that  guarantee  the  existence  of  a 
unique  polynomial  for  all  pennissible  combinations  of  n  and  m.  These  invariance 
properties  specify  that  the  polynomial  has  the  same  form  when  the  axes  are  rotated  about 
the  origin  by  a  continuous  function,  F,  with  a  period  of  2  n  and  an  angle  of  rotation  of  (f> 
[22] 

Z  (x,y)  =  F  (</>)Z  (x',yr)  (4.7) 

where 


x  =xcos^  +  ysin^ 
y' '  =  -x  sin  (f>  +  y  cos  (j) 


(4.8) 


It  is  important  to  note  that  as  stated  earlier,  Zernike  polynomials  are  normal  about 
a  continuous  unit  circle.  The  aperture  for  the  segmented  space  telescope  shown  in  Figure 
5,  or  any  Cassegrain  configuration  telescope  for  that  matter,  is  not  continuous  due  to  the 
hole  in  the  middle  of  the  primary  mirror  where  the  light  reflected  from  the  secondary 
mirror  passes  through.  In  order  to  account  for  this  discontinuity,  a  set  of  Zernike 
polynomials  has  been  modified  to  be  more  appropriate  for  use  on  a  unit  ring  instead  of  a 
unit  circle  [26].  For  this  research,  however,  the  fringe  Zernike  polynomials  were  used 
since  the  primary  goal  was  to  develop  a  general  technique.  Modification  to  the  annular 
polynomials  as  set  forth  in  [26]  is  left  for  future  work. 


For  lower  order  terms,  the  aberrations  described  by  the  Zernike  polynomials  take 
on  other  familiar  names,  for  example,  focus,  coma,  astigmatism,  and  so  on.  As  with  any 
series  expansion,  the  more  terms  used,  the  more  accurate  the  approximation  will  be.  For 
this  research,  2 1  Zernike  polynomials  are  used  to  represent  the  wavefront.  The  selection 
of  21  Zernikes  was  based  on  an  analysis  of  the  singular  values  of  the  original  system.  For 
comparison  purposes,  [12]  uses  8  Zernike  modes  for  wavefront  estimation  and  [13]  uses 
15  Zemikes.  Additional  research  could  be  performed  using  different  numbers  of 
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polynomials,  however,  that  effort  was  outside  the  scope  of  this  research.  Table  2 
contains  a  listing  of  these  21  polynomials  in  Cartesian  form.  For  completeness,  the 
piston  term  is  included,  but  it  is  a  non-controllable  mode  in  this  application.  For  those 
polynomials  associated  with  a  more  common  name,  the  common  name  is  also  given. 


# 

Polynomial 

Term 

0 

1 

Piston 

1 

X 

X-tilt 

2 

y 

Y-tilt 

3 

-1+2(x2+y2) 

Focus 

4 

2  2 

x  -y 

Astigmatism  plus  defocus 

5 

2xy 

Astigmatism  plus  defocus 

6 

-2x+3x(x2+y2) 

Coma  plus  tilt 

7 

-2y+3y(x2+y2) 

Coma  plus  tilt 

8 

1  -6(x2+y2)+6(x2+y2)2 

3rd  order  spherical  plus  focus 

9 

x3-3xy2 

5th  order  aberration 

10 

3x2y-y3 

5th  order  aberration 

11 

-3x2+3y2+4x2(x2+y2)-4y2(x2+y2) 

5th  order  aberration 

12 

-6xy+8xy(x2+y2) 

5th  order  aberration 

13 

3x-12x(x2+y2)+1 0x(x2+y2)2 

5th  order  aberration 

14 

3y-12y(x2+y2)+10y(x2+y2)2 

5th  order  aberration 

15 

1  +  12(x2+y2)-30(x2+y2)2+20(x2+y2)3 

5th  order  aberration 

16 

x4-6x2y2+y4 

7th  order  aberration 

17 

4x3y-4xy3 

7th  order  aberration 

18 

-4x3+1  2xy2+5x3(x2+y2)-1 5xy2(x2+y2) 

7th  order  aberration 

19 

-12x2y+4y3+15x2y(x2+y2)-5y3(x2+y2) 

7th  order  aberration 

20 

6x2-6y2-20x2(x2+y2)+20y2(x2+y2)+15x2(x2+y2)-1 5y2(x2+y2)2 

7th  order  aberration 

21 

1 2xy-40xy(x2+y2)+30xy(x2+y2)2 

7th  order  aberration 

Table  2  Table  of  Zernike  Polynomials  From  [24]. 

The  innovative  output  reduction  technique  using  Zemike  polynomials  is  based  on 
the  fact  that  the  majority  of  the  outputs  from  the  segmented  space  telescope  model  are 
slope  measurements  from  a  Shack-Hartmann  wavefront  sensor.  Based  on  the  work  of 
Southwell  in  [27],  these  slopes  can  be  related  to  the  partial  derivatives  of  Zernike 
polynomials  given  by 
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5 


(4.9) 


x 


I 


ak 


dzt(x>y) 

dx 


M 

v=2>, 


k= 1 


dzk(x>y) 

8y 


(4.10) 


where  Sx  and  Sy  are  the  x  and  y  slope  measurements,  Zk  are  the  Zernike  polynomials 


from  Table  2,  ak  are  the  expansion  coefficients  in  terms  of  the  Zernike  polynomials,  and 


M  is  the  number  of  polynomials  used,  in  this  case  2 1 .  The  Zernike  polynomials  are  a 
function  of  x  and  v,  so  they  must  be  evaluated  at  the  actual  lenslet  locations.  The  process 
for  doing  this  will  be  described  in  Chapter  V. 


Putting  Equations  (4.9)  and  (4.10)  into  matrix  form  yields 

S  =  dZa  (4.11) 

Pre-multiplying  both  sides  of  Equation  (4.11)  by  the  pseudo-inverse  of  the  matrix  of 
Zernike  partial  derivatives  yields  the  polynomial  coefficient  vector 

a  =  dZfS  (4.12) 

This  technique  is  used  for  the  first  time  in  this  research  as  a  means  of  performing 
model  reduction  to  facilitate  the  design  of  a  robust  controller.  Chapter  V  will  show  how 
this  was  accomplished  on  the  analytical  model.  The  original  system  output  of  Shack- 
Hartmann  slopes  has  been  replaced  by  an  output  of  Zernike  polynomial  coefficients.  In 
the  case  of  the  segmented  space  telescope  model,  what  was  originally  732  slope 
measurements  is  reduced  to  21  coefficient  outputs.  Clearly,  the  amount  of  reduction  is 
equal  to  the  number  of  Zernike  polynomials  used. 

As  will  be  discussed  in  Chapter  V,  the  SVD  method  can  be  used  to  reduce  both 
the  inputs  and  the  outputs  of  the  model.  In  this  research,  the  Zernike  polynomials  were 
only  used  to  reduce  the  number  of  outputs.  After  the  output  reduction  with  Zernike 
polynomials,  an  SVD  method  was  used  to  reduce  the  number  of  inputs.  Neither  of  these 
methods  has  any  affect  on  the  number  of  states  of  the  system,  which  is  addressed  next. 
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B. 


STATE  REDUCTION  WITH  HANKEL  SINGULAR  VALUES 


The  Hankel  singular  value  method  of  model  reduction  uses  a  balanced  stochastic 
approach,  described  in  [28],  [29],  [16],  and  [30].  It  begins  by  assuming  that  the  system 
model,  G,  is  stable  and  can  be  partitioned  as  shown 


G  = 


A 

B 

C 

D 

(4.13) 


From  this  model,  the  controllability  and  observability  grammians,  P  and  Q,  are 
computed  by  the  following  equations 


P  =  j  eA'BBTeAT,dt 


which  satisfies  the  Lyapunov  equation 

AP  +  PAt  +BBt  =0 
and 


satisfying 


Q  =  j  eAT,CTCeA'dt 


AtQ  +  QA  +  CtC  =  0 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


where  B  and  C  come  from  the  state  space  model.  The  Hankel  singular  values  are  then  the 
square  root  of  the  product  of  the  largest  eigenvalues  and  the  controllability  and 
observability  grammians  as  shown  by 

cth=4WQ)  (4.18) 

Next,  a  balanced  realization  is  created  to  balance  the  controllability  and 
observability  grammians.  A  nonsingular  transformation  matrix,  T,  is  introduced  such  that 


A 

B 

C 

D 

G  = 


which  results  in  the  following  controllability  and  observability  grammians 


TAT _1 

TB 

CT' 

D 

(4.19) 


P  =  TPT 


(4.20) 
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(4.21) 


q={t-1)t  qt~\ 

The  transformation  matrix,  T,  can  be  chosen  such  that  P  and  Q  are  diagonal  matrices 
with  identical  values  as  shown  by 

P  =  TPTt  =  X  (4.22) 

and 

Q  =  (t~1)T  QT~l  =  X  (4.23) 

where  X  are  the  Hankel  singular  values.  To  determine  T,  the  singular  value 
decomposition  is  used. 


The  Hankel  singular  values  can  be  placed  into  a  diagonal  matrix  such  that 


cr 


(*■).-  ^max^) 


0 


0  X 


(4.24) 


2  J 


These  Hankel  singular  values  can  then  be  plotted  to  determine  how  many  should  be 
retained  and  how  many  should  be  truncated.  The  states  that  are  truncated  are  the  least 
controllable  and  least  observables  states  of  the  system,  so  truncating  them  will  have 
minimal  effect  on  the  system. 


The  reduced  model,  GR  ,  can  be  shown  to  be  bounded  by  the  size  of  the  truncated 
singular  values  through  the  following  inequality 

\\0,-G\l<2tr(i:,)  (4.25) 

The  total  number  of  system  states  is  therefore  reduced  to  the  number  of  Hankel  singular 
values  retained. 
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V.  APPLICATION  TO  SEGMENTED  SPACE  TELESCOPE 

MODEL 


A.  MODEL  REDUCTION 

As  described  in  Chapter  II,  the  segmented  space  telescope  model  is  extremely 
large  and  complex.  Before  a  robust  controller  can  be  synthesized  for  the  model,  the 
number  of  inputs  and  outputs  must  be  reduced  to  ease  the  computational  burden.  The 
following  discussion  will  explain  two  different  model  reduction  techniques  used  to 
accomplish  this.  The  first  is  a  Singular  Value  Decomposition  reduction  for  both  inputs 
and  outputs.  The  second  technique  first  performs  a  Zernike  polynomial  reduction  on  the 
outputs  followed  by  a  Singular  Value  Decomposition  reduction  on  the  inputs.  The 
margins  and  performance  of  both  controllers  are  compared  against  each  other. 

1.  SVD  Only  Reduction 

Recall  that  the  space  telescope  model  is  a  state  space  model.  It  is  essentially 
comprised  of  four  separate  matrices,  which  define  a  standard  state  space  model  as  shown 
by 

x  =  Ax  +  Bit 

(5.1) 

y  =  Cx  +  Du 

A  Singular  Value  Decomposition  cannot  be  performed  on  a  state  space  model  but 
must  be  perfonned  on  an  individual  matrix.  Therefore,  a  single  matrix  representation  of 
the  output  is  defined  as 

y  =  H(s)u  (5.2) 

where 

H(s)  =  C(sl  -  A)1  B  +  D  (5.3) 

By  letting  s  go  to  zero,  the  steady  state  value  of  this  transfer  function  is  obtained  which 
neglects  any  dynamics  of  the  system 

H(0)  =  -CAlB  +  D  (5.4) 
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Since  no  dynamics  are  present,  this  transfer  function  can  also  be  thought  of  as  a 
reconstruction  matrix.  As  discussed  in  Chapter  IV,  the  SVD  of  H  yields  U,  £ ,  and  V 
matrices.  A  plot  of  the  singular  values  from  H  is  shown  in  Figure  22. 


Singular  Vaules  of  Original  System  Transfer  Function  Matrix 


Figure  22  Singular  Values  from  Original  System. 

The  number  of  rows  in  the  U  matrix  corresponds  to  the  total  number  of  sensor 
outputs  in  the  system.  The  column  vectors  represent  coefficients  of  the  entire  output 
vector.  To  reduce  the  number  of  system  outputs,  the  columns  of  U  can  be  truncated  to 
the  desired  number  of  outputs.  In  this  research,  2 1  outputs  were  desired,  so  only  the  first 
21  columns  of  U  were  retained.  The  reduced  size  of  U  is  therefore  [936  x  21].  Now, 
define  a  matrix  P  to  be  the  transpose  of  the  reduced  U  matrix 

P  =  [Ux..Uj  (5.5) 

At  this  point  a  reduced  system  output  vector,  y  ,  can  be  defined  as 

y  =  Py  (5.6) 

In  the  same  way  that  the  U  matrix  from  the  SVD  of  H  was  used  to  reduce  the 
number  of  system  outputs,  the  V  matrix  from  the  SVD  of  H  can  be  used  to  reduce  the 
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system  inputs.  V  has  997  rows,  corresponding  to  the  system  inputs.  By  again  selecting 
the  first  21  columns  of  V  and  truncating  the  rest,  we  can  define  a  matrix  M  with 
dimensions  [997  x  21]  as 

M=[Vv..V2l\  (5.7) 

The  total  system  input  vector,  u,  can  therefore  be  defined  as 

u  =  Mu  (5.8) 

where  u  is  the  reduced  input  vector.  Substituting  u  and  y  into  the  original  state  space 
model  in  Equation  (5.1)  yields  the  following  input/output  reduced  state  space  model 

x  =  Ax  +  BMu 

(5.9) 

y  =  PCx  +  PDMu 

This  model  is  rewritten  with  new  B  ,  C ,  and  D  matrices 


x  =  Ax  +  Bu 
y  =  Cx  +  Du 

where 


(5.10) 


B  =  BM 

C  =  PC  (5.11) 

D  =  PDM 

At  this  point,  the  number  of  inputs  and  outputs  in  the  system  have  been 
sufficiently  reduced  such  that  a  robust  controller  can  be  synthesized.  The  only  thing  that 
remains  is  to  reduce  the  number  of  states  of  the  system.  The  segmented  space  telescope 
model  has  332  states,  ft  is  possible  to  generate  a  robust  controller  based  on  a  system  with 
this  many  states;  however,  decreasing  the  number  of  states  will  facilitate  real  time 
implementation.  For  this  research,  a  state  reduction  was  performed  and  the  procedure  for 
that  reduction  will  be  discussed  after  first  applying  the  innovative  method  for 
input/output  reduction  using  Zemike  polynomials. 
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2. 


Zernike  Reduction 


The  process  of  performing  an  input/output  reduction  with  Zemike  polynomials  is 
actually  quite  similar  to  the  SVD  reduction  process  previously  discussed;  however,  to 
avoid  confusion  and  to  highlight  the  subtle  differences  the  entire  process  will  be 
presented.  As  before,  the  original  system  is  given  in  the  state  space  model  shown  in 
Equation  (5.1).  The  key  difference  in  this  reduction  method  occurs  at  this  point.  Before 
anything  else  is  done,  the  output  from  the  model  is  reduced  into  Zemike  polynomial 
coefficients  as  described  in  Chapter  IV. 

The  output  slopes  from  the  Shack-Hartmann  wavefront  sensor  are  equated  to  the 
partial  derivatives  of  the  Zemike  polynomials.  Any  number  of  Zernike  polynomials  may 
be  used;  2 1  were  used  for  this  research.  The  matrix  of  Zemike  partial  derivatives,  dZ,  is 
formed  with  dimensions  of  [732  x  21]  where  the  number  of  rows  corresponds  to  the 
number  of  Shack-Hartmann  slope  measurements  and  21  is  the  number  of  Zemike 
polynomials  used.  In  order  to  populate  the  dZ  matrix  with  actual  values,  the  physical 
locations  of  the  lenslets  in  the  model  are  used.  Refer  to  Figure  23,  which  shows  the 
lenslet  locations  and  numbering  scheme. 
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Figure  23  Mirror  Segment  and  Lenslet  Orientation. 


First,  assuming  the  physical  diameter  of  each  lenslet  to  be  1,  the  x  y  coordinates 
for  the  center  of  each  lenslet  in  segment  1  are  computed  in  relation  to  the  center  of  lenslet 
5,  which  is  at  the  center  of  the  segment.  For  example,  lenslet  4  has  coordinates  of 


(-1  0)  and  lenslet  15  has  coordinates  of 


1.5 


y/3y 


The  y  coordinates  for  the 


lenslets  in  this  segment  are  then  translated  a  distance  of  9  * 


73/ 


from  the  center  of  the 


entire  array,  O.  Next,  the  coordinates  for  the  remaining  five  segments  are  computed  by 

rotating  about  a  circle,  each  segment  being  displaced  -60°  from  the  previous  in  the 

segment  numbering  scheme.  Finally,  all  coordinates  are  scaled  to  the  actual  lenslet 

diameter  and  normalized  about  a  unit  circle.  Then  they  are  substituted  into  the  Cartesian 
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Zernike  polynomial  partial  derivative  equations,  which  populate  the  dZ  matrix.  The 
nonnalized  distance  from  the  center  of  the  array  is  the  same  as  that  for  lenslets  52  and  56. 

The  pseudo-inverse  of  the  matrix  dZ  is  then  used  to  define  a  new  output  vector  y 

y'  =  dZ'y  (5.12) 

The  dZ  matrix  formed  from  the  x  and  y  coordinates  of  the  analytical  model  has  a 
favorable  condition  number  of  9.98.  Recall,  however,  from  Chapter  II,  that  the  original 
system  model  contains  more  outputs  than  the  732  Shack-Hartmann  slopes.  There  are  936 
total  outputs,  which  include  measurements  from  other  sensors  as  well  as  the  system 
states.  Therefore,  a  slight  modification  is  made  to  the  pseudo-inverse  of  the  dZ  matrix  by 
adding  a  block  of  zeros  at  the  end  so  that  the  matrix  dimensions  match  those  of  the 
original  system.  The  modified  dZ '  is  shown  by 

Zernike  reduction  matrix  =  (  dZ'  |  0)  (5.13) 

where  the  dimensions  of  the  dZ'  matrix  are  [21  x  732]  and  the  block  of  zeros  has 
dimensions  of  [21  x  204]  yielding  a  final  dimension  of  [21  x  936]  for  the  Zernike 
reduction  matrix. 


At  this  point  the  output  reduced  system  model  has  the  following  form 

x  =  Ax  +  Bu 
y  =  dZ '  Cx  +  dZ '  Du 

which  can  be  rewritten  as 

x  =  Ax  +  Bu 
y  =  Cx  +  D’u 

where 


(5.14) 


(5.15) 


C  =  dZ^C  (5.16) 

D’  =  dZ^D  (5.17) 

Now  the  process  for  reducing  the  system  inputs  follows  very  closely  to  that  used  in  the 
SVD  reduction. 


A  transfer  function  matrix,  H,  is  now  defined  as 

H{s)  =  C\sl  -  A)1  B  +  D’ 


(5.18) 


56 


Letting  s  — >•  0 


H(0)  =  -C'A'B  +  D'  (5.19) 

The  SVD  for  this  transfer  function  of  the  output  reduced  system  model  is  taken, 
which  yields  the  U,  2 ,  and  V  matrices  as  before.  Figure  24  shows  the  singular  values 
from  this  matrix. 


Singular  Vaules  ofZernike  Output  Reduced  System  Transfer  Function  Matrix 


Figure  24  Singular  Values  from  Zernike  Output  Reduced  System. 

Since  the  system  has  already  been  output  reduced,  the  U  matrix  is  not  needed  this 
time,  and  only  the  V  matrix  is  used  to  reduce  the  inputs.  As  before,  21  columns  of  V  are 
retained  so  that  Fhas  dimensions  of  [997  x  21].  Defining  M'  to  be 

M'  =  [V ;...F21]  (5.20) 

the  total  system  input  vector,  u,  can  therefore  be  defined  as 

u  =  M'u' 


Substituting  this  into  the  Zernike  reduced  system  model  yields 


(5.21) 


(5.23) 


This  model  can  be  simplified  as  before  to 

x  =  Ax  +  B'u 
y  =  C'x  +  D'u' 

where 

B'  =  BM'  (5.24) 

and  D'  is  modified  from  Equation  (5.17)  to  be 

D'  =  dZfDM'  (5.25) 

At  this  point,  the  original  system  state  space  model  has  had  its  inputs  and  outputs 
reduced  by  two  different  methods:  an  SVD  only  method  to  reduce  both  inputs  and 
outputs;  and  a  Zernike  polynomial  method  to  reduce  the  outputs,  followed  by  an  SVD 
method  to  reduce  the  inputs.  From  this  point,  the  remaining  steps  in  designing  a  robust 
controller  are  identical. 

3.  State  Reduction 

The  original  system  state  space  model  contains  332  states.  Neither  of  the 
input/output  reduction  techniques  discussed  has  changed  or  altered  the  number  of  states. 
As  mentioned  above,  a  state  reduction  is  not  needed  to  generate  a  controller;  however,  the 
designed  controller  will  have  the  same  number  of  states  as  the  system  it  is  designed  to 
control.  Therefore,  it  may  be  advantageous  to  reduce  the  number  of  system  states  prior  to 
generating  the  controller  in  order  to  facilitate  real  time  implementation  on  a  system. 

Equations  (5.10)  and  (5.23)  both  represent  input/output  reduced  state  space 
models,  one  for  each  reduction  technique  used.  From  this  point  on,  the  discussion  will 
use  the  nomenclature  used  in  Equation  (5.10),  although  nomenclature  from  either 
equation  would  be  suitable. 

The  number  of  system  states  is  reduced  using  the  Hankel  Singular  Value  state 
reduction  techniques  discussed  in  Chapter  IV.  A  plot  of  the  Hankel  Singular  Values  from 
the  input/output  reduced  system  model  is  shown  in  Figure  25. 
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Figure  25  Flankel  Singular  Values. 


Based  on  this  graph,  the  decision  was  made  to  reduce  the  number  of  system  states 
to  240,  thereby  retaining  the  most  significant  singular  values.  By  defining  x  as  the 
Flankel  reduced  state,  the  new  system  state  space  model  can  be  written  as 


(5.26) 


where  A',  B' ,  C',  and  y'  all  reflect  the  reduced  number  of  states.  This  system  contains 
240  states,  2 1  inputs,  and  2 1  outputs.  From  here  it  is  possible  to  proceed  with  the  robust 
controller  synthesis. 

B.  ROBUST  CONTROLLER  SYNTHESIS 

After  the  model  is  sufficiently  reduced  to  allow  synthesis  of  a  robust  controller, 
the  next  step  in  the  process  is  to  generate  the  2-input,  2-output  construct  as  shown  in 
Figure  19.  Part  of  this  process  involves  incorporating  the  weights,  Wi  and  W$.  Wj  is  the 
weight  on  the  error  between  the  reference  signal  and  measured  output;  Wj  is  the  weight 
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on  the  measured  output.  For  this  design,  first  order  transfer  functions  were  chosen  for 
each  and  are  shown  in  Equation  (5.27)  and  Equation  (5.28).  Figure  26  shows  a  Bode  plot 
of  Wi  and  the  inverse  of  W3. 


0.55  +  17.32 
5  +  0.3465 


(5.27) 


105  +  2095 
5  +  2116 


Weighting  funtion  Bobe  plots 


Figure  26  Model  Weighting  Function  Bode  Plot. 


Each  weighting  function  specifies  a  low  frequency  gain,  a  crossover  frequency, 
and  a  high  frequency  gain.  The  effect  is  analogous  to  classical  loopshaping  techniques. 
These  transfer  functions  were  determined  based  on  their  ability  to  achieve  the  desired 
bandwidth  for  the  closed  loop  system,  while  at  the  same  time  not  being  so  restrictive  as  to 
preclude  the  synthesis  of  a  robust  controller.  The  desired  control  bandwidth  for  the  space 
telescope  model  was  10  Hz.  The  classical  controller  designed  for  the  model  was  unable 
to  achieve  this  bandwidth,  which  was  a  driving  motivation  for  this  research  into  robust 
control  techniques. 
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The  result  of  augmenting  the  reduced  system  model  with  these  weighting 
functions  is  a  system  in  the  form  shown  in  Figure  19.  The  total  number  in  inputs  in  this 
augmented  system  is  42,  resulting  from  21  inputs  from  the  reference  signal  and  21  inputs 
from  the  control  signal.  There  are  63  total  outputs  in  the  augmented  system  resulting 
from  21  error  signal  outputs  multiplied  by  Wj,  21  measured  outputs  multiplied  by  Wj,  and 
21  error  signal  outputs  sent  to  the  controller.  The  augmented  model  has  440  states.  The 
Hx  robust  controller  that  is  synthesized  for  this  system  has  21  inputs,  21  outputs,  and 

440  states  regardless  of  whether  the  model  reduction  is  performed  with  the  SVD  only 
method  or  the  Zernike  polynomial  /  SVD  combination.  This  controller  can  now  be  used 
to  check  the  performance  on  the  original  plant  model.  The  results  from  these  simulations 
are  given  in  Chapter  VI. 

C.  PERFORMANCE  AND  STABILITY  MARGINS 

The  first  analysis  of  the  quality  of  the  two  robust  controllers  designed  is 
accomplished  by  examining  the  poles  of  the  controllers.  Figure  27  and  Figure  28  show 
the  location  of  the  poles  for  the  closed  loop  system  from  each  controller.  The  closed  loop 
system  combines  the  Hx  controller  with  the  full  plant  and  assumes  negative  unity 
feedback.  The  poles  for  each  controller  are  plotted  on  the  same  scale. 
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x  -]g4  Closed  loop  poles  for  SVD  only  reduction 


Figure  27  Closed  Loop  Poles  for  SVD  only  Reduced  Model. 


x  i  q4  Closed  loop  poles  for  Zernike  reduction 
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Figure  28  Closed  Loop  Poles  for  Zernike  Polynomial  Reduced  Model. 
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The  following  show  the  same  results  zoomed  in  on  the  area  near  the  imaginary  axis. 


Imaginary  axis  vicinity  for  SVD  reduction  closed  loop  poles 


Figure  29  Imaginary  Axis  Vicinity  for  SVD  Reduction  Closed  Loop  Poles. 


Imaginary  axis  vicinity  forZernike  reduction  closed  loop  poles 


Figure  30  Imaginary  Axis  Vicinity  for  Zernike  Polynomial  Reduction  Closed  Loop 

Poles. 
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As  can  be  seen  in  the  plots,  the  controller  for  the  Zernike  polynomial  reduced 
model  has  poles  much  closer  to  the  imaginary  axis;  however,  in  both  cases  all  poles  are  in 
the  left  half-plane,  indicating  the  stability  of  the  closed  loop  systems.  Table  3 
summarizes  the  poles  that  are  closest  to  the  imaginary  axis  for  each  of  the  model 
reduction  techniques  for  both  closed  loops  systems. 


SVD  Reduced  Model 

-0.7118  ±268.69/ 
-0.8729  ±4545.6023/ 
-0.9149  ±  182.9182/ 


Zemike  Reduced  Model 

-0.0327  ±406.1887/ 
-0.0623  ±  1098.7036/ 
-0.6178  ±  588.1804/ 


Table  3  Closed  Loop  Poles  Closest  to  jco  Axis. 

The  next  level  of  analysis  looks  at  the  singular  values  of  the  plant,  the  controllers, 
and  the  closed  loop  systems.  Recall  that  a  singular  value  plot  is  the  MIMO  equivalent  to 
the  Bode  magnitude  plot  of  a  SISO  system.  Figure  31  shows  the  singular  values  for  the 
input/output  reduced  plant  with  all  332  states.  Computing  the  singular  values  on  the  full 
plant  with  all  inputs  and  outputs  was  too  computationally  intensive. 


Input/Output  Reduced,  332  State  Plant  Singular  Values 


Figure  3 1  Singular  Values  for  Original  System  Model,  Gred  . 
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As  evident  in  this  plot,  the  main  natural  frequencies  of  the  plant  occur  between  3 1 
and  760  Hz.  A  properly  designed  controller  will  provide  stability  and  negate  any  adverse 
effects  due  to  these  resonances  in  achieving  the  desired  closed  loop  performance. 

The  singular  values  for  the  controller  from  the  SVD  only  reduced  model  and  from 
the  Zernike  polynomial  reduced  model  are  shown  in  Figure  32  and  Figure  33.  The 
controller  is  represented  by  KSVD  or  K/ern  depending  on  how  the  model  reduction  was 

performed.  In  both  cases,  the  plots  have  been  zoomed  in  on  the  applicable  range  of 
frequencies  between  31  and  760  Hz  where  the  controllers  have  been  designed  to  handle 
the  resonances.  Based  on  a  visual  comparison,  it  appears  that  both  controllers  should  be 
able  to  compensate  for  the  resonances  in  the  system  and  provide  stability.  This 
observation  will  be  verified  for  the  closed  loop  system. 


SVD  Reduced  Controller  Singular  Values 


Frequency  (Hz) 


Figure  32  Singular  Values  for  H  a  Controller,  KSVD . 
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Zernike  Reduced  Controller  Singular  Values 


Figure  33  Singular  Values  for  Hx  Controller,  KZem . 

The  next  analysis  is  to  examine  the  open  loop  singular  values  shown  in  Figure  34 
and  Figure  35.  The  open  loop  system  is  calculated  as  either  GKSVD  or  GKZern .  Notice  in 

these  two  plots  that  the  curve  is  smoother  for  the  controller  based  on  the  Zernike  reduced 
model.  There  are  fewer  residual  peaks  left  from  the  natural  frequencies  of  the  system. 
This  is  an  indication  that  the  controller  for  the  Zernike  reduced  model  more  closely 
matches  the  original  system.  It  is  better  able  to  compensate  for  the  plant’s  natural 
frequencies.  While  not  an  indication  of  performance  or  robustness,  this  is  an  indication 
that  the  controller  should  be  stable. 
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Singular  Values  (dB)  «  Singular  Values  (dB' 


SVD  Reduced  Open  Loop  System  Singular  Values 


rigure  34  Open  Loop  System  Singular  Values  for  GK 


SVD  ' 


Zernike  Reduced  Open  Loop  System  Singular  Values 


Figure  35  Open  Loop  System  Singular  Values  for  GKZern . 
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Next,  the  closed  loop  system  is  formed  by  providing  a  unity  negative  feedback 
loop  to  the  open  loop  systems  GKSVD  and  GKZem  as  shown  by 


y  = 


GK 


SVD 


\  +  GK 


SVD 


(5.29) 


y  = 


GK 


Zern 


1  +  gk7 


(5.30) 


The  singular  value  plots  from  both  closed  loop  systems  are  shown  in  Figure  36  and 
Figure  37. 


SVD  Reduced  Closed  Loop  System  Singular  Values 
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0 


Zernike  Reduced  Closed  Loop  System  Singular  Values 


As  with  the  open  loop  singular  value  plots,  the  smoother  shape  of  the  Zemike 
reduced  controller  closed  loop  singular  value  plot  indicates  that  it  more  closely  matches 
the  natural  frequencies  of  the  original  plant.  This  observation  will  be  checked  in  the 
simulation  analyses. 

The  shape  of  the  open  and  closed  loop  singular  value  plots  for  both  controllers 
matches  the  desired  shape  specified  by  the  weighting  functions  shown  in  Figure  26.  The 
weighting  functions  specified  a  low  frequency  performance  requirement  with  a  gain  no 
larger  than  5  dB  at  10  FIz.  It  also  satisfies  the  high  frequency  robustness  requirement  of 
no  more  than  -15  dB  at  1000  Hz. 

Also  evident  in  the  singular  value  plots  is  the  control  bandwidth  for  each 
controller.  The  bandwidth  was  determined  by  identifying  where  the  singular  value  plots 
crossed  the  -3dB  line  for  the  closed  loop  system.  The  SVD  reduced  model  controller  has 
a  control  bandwidth  in  the  range  of  20-24  Hz,  whereas  the  Zernike  reduced  model 
controller  bandwidth  resides  from  15-19  Hz.  Both  of  these  bandwidths  are  higher  than 
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that  obtained  using  classical  control  methods  leading  to  the  conclusion  that  the  robust 
control  techniques  employed  in  the  design  process  provide  better  performance  and  are 
better  able  to  accomplish  disturbance  rejection. 

D.  ANALYTICAL  RESULTS 

This  chapter  will  present  simulation  results  for  the  two  different  robust  controllers 
developed  for  the  segmented  telescope  model.  One  controller  was  designed  based  on  a 
model  whose  inputs  and  outputs  were  both  reduced  by  the  SVD  technique;  the  model  for 
the  other  controller  had  its  outputs  reduced  with  Zemike  polynomials  and  inputs  reduced 
with  the  SVD  technique.  The  Hankel  singular  value  state  reduction  and  all  other  aspects 
of  controller  design  were  identical. 

A  simple  Simulink  model  shown  in  Figure  38  was  used  to  validate  the 
performance  of  both  controllers. 


Figure  38  Segmented  Space  Telescope  Controller  Validation  Model. 

The  full  plant  used  for  both  cases  was  the  original  997  input,  936  output,  332  state 
model.  The  controller  for  both  cases  had  21  inputs  and  21  outputs.  The  number  of  states 
in  the  controller  does  not  have  to  match  the  number  of  states  in  the  original  plant  as  long 
as  the  inputs  and  outputs  go  through  the  appropriate  projection.  The  full  output  that  was 
sent  to  the  workspace  for  plotting  consisted  of  only  Shack-Hartmann  slope  data.  The 
phase  diversity  and  other  outputs  were  discarded  and  not  plotted.  The  actuator  projection 
is  the  SVD  reduction  of  the  inputs,  which  consists  of  the  first  2 1  columns  of  the  V  matrix 
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obtained  from  the  SVD.  This  recreates  the  full  control,  u,  from  the  modified  control,  u  , 
generated  from  the  reduced  controller.  Equation  (5.8)  describes  this  process. 

For  the  case  of  the  Zernike  polynomial  reduced  model,  the  following  conditions 
apply.  The  reference  input  vector  and  the  reduced  output  sent  to  the  workspace  were 
both  the  Zernike  polynomial  coefficients.  The  sensor  projection  is  the  pseudo-inverse  of 
the  dZ  matrix  formed  from  the  partial  derivatives  of  the  Zernike  polynomials.  The 
reduced  output  becomes  that  as  shown  in  Equation  (5.12). 

For  the  case  of  the  SVD  only  reduced  model,  the  following  conditions  apply.  The 
reference  input  vector  and  the  reduced  output  sent  to  the  workspace  were  both  the 
coefficients  of  the  full  output  vector.  The  sensor  projection  is  the  transpose  of  the  first  21 
columns  of  the  U  matrix  formed  from  the  SVD  of  the  model.  The  reduced  output 
becomes  that  as  shown  in  Equation  (5.6). 

The  band-limited  white  noise  was  used  to  add  disturbances  to  the  simulation.  In 
some  cases  no  disturbance  was  added,  but  when  it  was  added,  it  provided  a  random  step 
to  all  control  inputs  at  a  specified  time  interval.  Figure  39  shows  what  this  would  look 
like.  This  type  of  disturbance  can  be  thought  of  as  a  changing  initial  condition  as  would 
be  encountered  if  the  space  telescope  performed  a  slewing  maneuver  to  acquire  a  new 
target. 
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Figure  39  Disturbance  Input  Representation. 

1.  System  Performance 

The  first  test  of  the  controller  performance  was  how  well  the  system  could  track  a 
step  input  of  0.2  on  all  reference  inputs.  Figure  40  and  Figure  41  show  the  results  for 
both  controllers. 


Disturbance  Input  Steps 
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Amplitude  Amplitude 


0.2  Step  Response  for  SVD  Reduced  Controller 


Figure  40  Step  Input  for  SVD  Only  Reduced  Model. 


0.2  Step  Response  for  Zernike  Reduced  Controller 


Figure  41  Step  Input  for  Zernike  Polynomial  Reduced  Model. 
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These  plots  were  only  used  to  verify  that  each  controller  is  capable  of  achieving  a 
desired  steady  state  response  to  a  specified  reference  input.  In  this  case,  both  controllers 
are  capable  of  accomplishing  this.  No  attempt  was  made  at  comparing  the  perfonnance 
of  the  two  controllers  based  on  these  plots,  however,  since  the  output  is  different  in  each 
case.  For  the  SVD  reduced  controller,  the  output  is  a  coefficient  of  a  total  output  vector, 
whereas  the  output  of  the  Zemike  reduced  controller  is  a  Zemike  polynomial  expansion 
coefficient. 

2.  Disturbance  Rejection 

Now  that  both  controllers  have  been  shown  to  be  stable  and  can  achieve  a  desired 
steady  state,  it  is  time  to  compare  their  performance.  The  perfonnance  criterion  is  to 
achieve  a  planar  wavefront,  or  in  other  words,  the  slopes  from  the  Shack-Hartmann 
wavefront  sensor  should  approach  zero.  An  equivalent  statement  is  that  the  expansion 
coefficients  of  the  Zernike  polynomials  all  become  zero.  If  the  Zemike  polynomial 
expansion  coefficients  are  zero  then  there  are  no  aberrations  in  the  wavefront  with  the 
possible  exception  of  a  piston  component,  which  is  not  measurable  in  this  setup.  All 
other  aberrations  should  be  zero. 

In  order  to  make  a  valid  comparison,  the  output  from  each  controller  needs  to 
represent  the  same  data.  Recall,  however,  that  the  output  from  the  SVD  reduced 
controller  was  not  a  Zernike  polynomial  expansion  coefficient,  but  rather  a  coefficient  on 
the  total  output.  This  difference  can  be  remedied  by  converting  the  full  slope  output  from 
the  SVD  reduced  controller  in  a  Zemike  polynomial  expansion  coefficient.  This  is  done 
by  multiplying  the  full  output  by  the  pseudo-inverse  of  the  dZ  matrix  from  Equation 
(5.12).  In  the  simulation  block  diagram,  this  is  represented  by  the  extra  block  in  green  in 
Figure  42. 
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Figure  42  Simulation  Model  with  Zernike  Coefficient  Conversion. 

The  simulation  was  run  for  4.1  seconds  with  a  reference  input  of  zero.  The  band 
limited  white  noise  disturbance  represented  in  Figure  39  was  set  to  provide  a  new  random 
step  disturbance  every  second.  The  output  from  the  SVD  reduced  controller,  converted 
to  Zemike  polynomial  expansion  coefficients,  is  shown  in  Figure  43. 


Zernike  Coefficients  from  SVD  Reduced  Controller  with  Random  Step  Disturbances 


Figure  43  Zemike  Polynomial  Coefficients  in  Presence  of  Random  Disturbance  for  SVD 

Reduced  Controller. 
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Notice  how  the  controller  does  a  good  job  of  bringing  the  disturbance  back  close 
to  zero  in  about  0.5  seconds.  However,  some  of  the  coefficients  asymptotically  go  to 
other  non-zero  values  and  exhibit  oscillatory  behavior  around  those  values.  In  some 
cases,  the  coefficients  are  off  by  as  much  as  ±  0.04.  This  result  can  be  compared  with 
that  of  the  Zernike  reduced  controller  shown  in  Figure  44. 

Zernike  Coefficients  from  Zernike  Reduced  Controller  with  Random  Step  Disturbances 


Figure  44  Zernike  Polynomial  Coefficients  in  Presence  of  Random  Disturbance  for 

Zernike  Reduced  Controller. 

Notice  here  that  in  roughly  the  same  amount  of  time  the  Zernike  reduced 
controller  drives  all  Zernike  polynomial  coefficients  to  a  zero  value  with  negligible 
deviation  from  zero  and  negligible  oscillatory  behavior  in  the  steady  state.  The 
conclusion  from  this  simulation  is  that  the  Hx  robust  controller  synthesized  for  the 
Zernike  polynomial  reduced  model  is  better  at  achieving  a  planar  wavefront  in  the 
presence  of  disturbances.  This  is  not  an  unexpected  conclusion,  since  the  Zernike 
polynomial  reduction  was  theorized  to  be  uniquely  suited  to  reduce  the  Shack-Hartmann 
wavefront  sensor  outputs  of  the  segmented  space  telescope  model. 
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VI.  EXPERIMENTAL  APPLICATION  TO  ADAPTIVE  OPTICS 

TESTBED 


A.  SYSTEM  IDENTIFICATION 

This  chapter  will  describe  the  work  to  adapt  the  robust  control  techniques  for  use 
on  the  Spacecraft  Research  and  Design  Center’s  adaptive  optics  testbed,  as  described  in 
Chapter  II.  As  noted  in  the  literature  [9],  this  is  not  the  first  use  of  a  robust  controller 
designed  for  a  laboratory  testbed.  The  intent  here  is  not  to  replicate  the  previous  work 
done,  but  to  validate  the  design  process  used  for  the  analytical  model  on  testbed 
hardware. 

To  do  this,  a  robust  controller  was  designed  for  the  first  deformable  mirror  / 
wavefront  sensor  closed  loop  system  on  the  SRDC  testbed.  It  was  originally  thought  that 
since  the  first  deformable  mirror  was  an  extremely  thin  membrane  mirror,  it  would  be 
lightweight  and  prone  to  have  low  natural  frequencies  just  as  a  lightweight  segmented 
space  telescope  mirror  would.  This  would  allow  verification  of  the  robust  control 
techniques  on  systems  with  unmodeled  system  dynamics. 

The  difficulty  is  that  no  model  of  the  actual  hardware  exists.  Therefore,  it  was 
necessary  to  perform  a  system  identification  on  the  mirror  to  detennine  its  natural 
frequencies.  Using  a  position  sensing  device  it  was  determined  that  the  lowest  natural 
frequency  was  around  918  Hz,  much  higher  than  originally  thought.  It  is  believed  that 
the  electrostatic  actuators  provide  a  degree  of  rigidity  to  the  system,  thus  reducing  the 
amount  of  dynamics  that  might  be  present  in  the  system.  In  the  design  of  actual  flight 
systems,  this  would  actually  be  a  benefit  as  there  would  be  less  uncertainty  in  the  system 
that  the  controller  must  be  able  to  handle.  For  the  purposes  of  this  research,  however,  it 
meant  that  a  new  approach  had  to  be  developed. 

B.  ROBUST  CONTROLLER  SYNTHESIS 

Since  the  natural  frequencies  of  the  membrane  deformable  mirror  were  around 
918  Hz,  it  was  decided  that  the  mirror  could  be  modeled  as  a  static  system  and  a  robust 
controller  designed  for  this  static  system.  During  real  time  wavefront  correction, 
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disturbance  inputs  could  be  injected  into  the  system  and  the  robust  controller 
performance  measured  in  the  presence  of  these  disturbances. 

The  best  choice  for  a  static  system  model  is  the  reconstruction  or  influence  matrix 
of  the  system.  Recall  that  the  reconstruction  matrix  is  obtained  by  actuating  (poking) 
each  of  the  mirror  actuators  in  turn,  allowing  sufficient  settling  time,  and  measuring  the 
resulting  effect  on  the  wavefront  with  the  Shack-Hartmann  wavefront  sensor.  In  classical 
control  applications,  the  wavefront,  (f> ,  is  measured  and  multiplied  by  the  inverse  (or  in 
this  case  the  pseudo-inverse)  of  the  reconstruction  matrix,  R ,  to  determine  the  control 
input,  u 

u  =  Rf</>  (6.1) 

In  the  case  of  the  membrane  defonnable  mirror  on  the  Spacecraft  Research  and 
Design  Center  AO  testbed,  the  control  input  is  a  vector  of  37  actuator  voltages,  one  for 
each  actuator.  The  wavefront  measurement  is  a  vector  of  254  slopes,  two  for  each  Shack- 
Hartmann  lenslet.  The  reconstruction  matrix  therefore  has  dimensions  of  [254  x  37]. 
The  large  difference  between  the  number  of  rows  and  columns  in  this  matrix  is  one  factor 
that  increases  the  condition  number  to  the  rather  high  value  of  33,628.  This  can  lead  to 
numerical  error  when  taking  the  pseudo-inverse  of  this  matrix  by  the  Moore-Penrose 
method  as  given  by 

Xf  =XT[XXT)~1  (6.2) 

What  is  needed  then  is  a  way  to  reduce  the  number  of  measurements.  The 
previously  discussed  method  of  model  reduction  with  Zernike  polynomials  is  quite 
appropriate  for  this  task.  By  using  the  Zernike  polynomial  output  reduction  method 
employed  on  the  analytical  model  and  projecting  the  254  slope  measurements  into  21 
Zernike  polynomials,  the  resulting  vector  has  only  21  coefficients.  This  reduced 
reconstruction  matrix  now  has  dimensions  of  [21  x  37],  with  a  reduced  condition  number 
of  257  and  a  less  complex  model  on  which  to  design  a  robust  controller.  This  reduced 
reconstruction  matrix  is  used  with  the  following  weighting  functions  to  generate  the 
robust  controller 

0.55  +  17.32 

Rj  = 


5  +  0.3465 
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(6.3) 


(6.4) 


fK  = 


10^  +  2095 


5  +  2116 

A  plot  of  these  weighting  functions  is  show  in  Figure  45. 


Weighting  funtion  Bobe  plots 


Figure  45  AO  Testbed  Robust  Controller  Weighting  Functions. 

The  resulting  Hx  controller  is  a  state  space  model  with  2 1  inputs,  which  are  the 

Zernike  polynomial  expansion  coefficients,  37  outputs,  which  are  translated  into  actuator 
voltages,  and  42  states.  When  the  controller  is  implemented  on  the  testbed,  these  42 
states  must  be  integrated  for  each  iteration  through  the  control  loop.  In  order  to  minimize 
the  computational  effort  associated  with  all  these  integrations,  the  continuous 
controller  is  discretized  for  implementation  according  to  the  following  set  of  equations 

x(&  + 1)  =  Ax{k)  +  Bu{k) 
y(k)  =  Cx(A:)  +  Di/(k) 

The  singular  value  plot  of  the  AO  testbed  controller  is  shown  in  Figure  46. 
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AO  Testbed  Controller  Singular  Values 


Figure  47  and  Figure  48  show  the  testbed  open  loop  and  closed  loop  singular  values. 


AO  Testbed  Open  Loop  Singular  Values 
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AO  Testbed  Closed  Loop  Singular  Values 


The  smoothness  of  these  plots  is  attributable  to  the  fact  that  there  are  no  dynamics 
or  resonances  in  the  model  used.  There  is  nothing  in  these  plots  to  indicate  that  the 

controller  designed  for  the  testbed  will  be  ineffective  at  correcting  the  wavefront 
aberrations.  The  experimental  results  will  show  this  to  be  the  case. 

C.  EXPERIMENTAL  RESULTS 
1.  System  Performance 

The  results  of  the  AO  testbed  controller  will  be  compared  with  the  results 
obtained  using  a  classical  reconstruction  matrix  approach.  The  wavefront  estimation  was 
performed  by  relating  the  Shack-Hartmann  wavefront  sensor  slope  measurements  to  the 
partial  derivatives  of  the  Zemike  polynomials  as  shown  in  Equations  (4.9)  and  (4.10). 
The  control  law  used  an  iterative  feedback  loop  based  on  the  reconstruction  matrix  as 
developed  in  [3]  and  given  as 

H,+i  ='hl~gR\  (6-6) 
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where  g  is  a  gain,  and  a  is  the  vector  of  Zernike  polynomial  expansion  coefficients  for 
each  iteration,  n.  Before  describing  the  performance,  it  is  necessary  to  describe  the 
evaluation  criteria  used. 


The  first  performance  measure  used  is  the  peak-to-valley  wavefront  aberration. 
This  figure  of  merit  is  a  sum  of  the  absolute  values  of  the  minimum  and  maximum 
deviation  from  a  planar  wavefront.  The  resulting  number  is  a  representation  of  the 
number  of  waves  of  aberrations  present.  The  piston,  tip,  and  tilt  aberrations  present  in  the 
wavefront  are  neglected  in  order  to  allow  examination  of  the  underlying  higher  order 
aberrations.  The  peak-to-valley  wavefront  error  is  then  represented  in  a  three- 
dimensional  circular  plot. 


The  second  perfonnance  measure  is  the  root  mean  square  (RMS)  error.  This 
figure  of  merit  is  the  square  root  of  the  wavefront  variance  and  is  given  as 

j  i  VF7  2 

(7l=  —  \  I  Mx’-f)-$>(x’-v)]  dxdy  (6.7) 

This  variance  is  related  to  the  Zernike  polynomial  coefficients  through  the 
following  relationship,  where  the  reference  wavefront,  </>0 ,  is  assumed  to  be  planar 

M 

(/){x,y)-(j)Q  (x,y)  =  Y,akZk  ix>y)  (6-8) 

k= 1 

The  RMS  error  can  then  be  written  as 


O’, 


M 


1  1  \ll-x2 

t!  1  \_Zk  (x,y)]  dxdy 


where  ak  has  been  pulled  in  front  of  the  integral  because  it  is  a  constant. 


(6.9) 


The  control  loop  was  run  for  100  iterations  with  each  controller.  At  the  end  of 
each  run,  the  wavefront  was  sampled  with  a  science  camera  to  determine  the  final 
wavefront  peak-to-valley  error.  A  representation  of  the  wavefront  for  each  simulation 
was  plotted.  The  RMS  error  for  each  iteration  was  also  plotted.  The  results  of  the 
classical  controller  are  shown  in  Figure  49  and  Figure  50.  A  gain  value  of  0. 1  was  used. 
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Wavefront  Surface,  Peak  to  Valley  =  4.7498 
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Figure  49  Classical  Controller  Peak-to-Valley  Wavefront  Error. 
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Figure  50  Classical  Controller  RMS  Wavefront  Error. 


RMS  Wavefront  Error  vs.  Iterations,  a  =  4.1024 
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As  can  be  seen,  a  steady  state  wavefront  error  was  obtained  after  about  25  iterations. 

Next,  the  robust  controller  was  used  with  all  the  same  settings  and  parameters. 
The  only  change  was  the  replacement  of  the  line  of  code  representing  the  classical 
controller  in  Equation  (6.6)  with  the  code  to  represent  the  discrete  robust  controller  in 
Equation  (6.5).  The  wavefront  estimation  algorithms  were  identical.  The  peak-to-valley 
and  RMS  results  from  100  iterations  with  the  robust  controller  are  shown  in  Figure  51 
and  Figure  52. 

Wavefront  Surface,  Peak  to  Valley  =  4.473 
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Figure  5 1  Robust  Controller  Peak-to- Valley  Wavefront  Error. 
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RMS  Wavefront  Error  vs.  Iterations,  a=  3.878 


Figure  52  Robust  Controller  RMS  Wavefront  Error. 


As  seen  in  the  plots,  the  robust  controller  outperforms  the  classical  controller  in 
peak-to-valley  wavefront  aberration,  final  RMS  value,  and  the  number  of  iterations 
required.  Increasing  the  gain  on  the  classical  controller  can  improve  the  response  time, 
but  this  also  leads  to  instability.  Other  attempts  to  design  robust  controllers  for  AO 
testbeds  have  yielded  lower  absolute  peak-to-valley  and  RMS  values  [9].  The  purpose 
here,  however,  was  not  to  achieve  the  lowest  absolute  error,  but  rather  to  experimentally 
verify  the  robust  control  synthesis  procedures  used  on  the  analytical  model.  The 
improved  performance  of  the  robust  controller  over  the  classical  controller  is  taken  as 
verification  of  this  approach.  The  higher  error  values  obtained  here  are  most  likely  the 
result  of  a  rudimentary  wavefront  estimation  algorithm.  Improving  the  centroiding 
algorithm  used  for  wavefront  estimation  is  left  for  future  work,  but  that  work  was  outside 
the  scope  of  this  research. 

2.  Disturbance  Rejection 

Additional  comparison  simulations  were  run  on  the  AO  testbed  with  a  disturbance 
on  the  control  voltages  applied  to  the  deformable  mirror.  The  disturbance  used  was  a 
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sinusoidal  input  at  varying  frequencies.  In  each  test  case,  the  robust  controller  proved 
superior  at  suppressing  the  disturbance  input  while  achieving  the  results  obtained  from 
the  simulation  with  no  sinusoidal  disturbance.  The  following  set  of  figures  show  results 
for  each  controller  for  two  different  disturbance  frequencies,  1.3963  Hz  and  6.2832  Hz. 
The  final  error  results  from  the  simulations  are  summarized  in  Table  4. 


Wavefront  Surface,  Peak  to  Valley  =  8.5044 


Figure  53  Classical  Controller  Peak-to- Valley  Wavefront  Error  with  1 .3963  Hz 

Disturbance. 
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RMS  Wavefront  Error  vs.  Iterations,  o=  4.6837 


Figure  54  Classical  Controller  RMS  Wavefront  Error  with  1.3963  Hz  Disturbance. 


Wavefront  Surface,  Peak  to  Valley  =  4.5052 
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Figure  55  Robust  Controller  Peak-to-Valley  Wavefront  Error  with  1.3963  Hz 

Disturbance. 
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RMS  Wavefront  Error  vs.  Iterations,  o  =  3.9232 


Figure  56  Robust  Controller  RMS  Wavefront  Error  with  1 .3963  Hz  Disturbance. 

Wavefront  Surface,  Peak  to  Valley  =  5.5255 
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Figure  57  Classical  Controller  Peak-to- Valley  Wavefront  Error  with  6.2832  Hz 

Disturbance. 
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RMS  Wavefront  Error  vs.  Iterations,  a  =  4.31 12 


Figure  58  Classical  Controller  RMS  Wavefront  Error  with  6.2832  Hz  Disturbance. 


Wavefront  Surface,  Peak  to  Valley  =  4.4815 
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Figure  59  Robust  Controller  Peak- to- Valley  Wavefront  Error  with  6.2832  Hz 

Disturbance. 
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RMS  Wavefront  Error  vs.  Iterations,  o  =  3.9349 


Figure  60  Robust  Controller  RMS  Wavefront  Error  with  6.2832  Hz  Disturbance. 


1.3963  Hz  Disturbance 

Peak-to-Vallev  Wavefront  Error 

RMS  Error  (o) 

Iterations 

Classical  Controller 

8.5044 

4.6837 

8 

Robust  Controller 

4.5052 

3.9232 

4 

6.2832  Hz  Disturbance 


Classical  Controller 

5.5255 

4.3112 

20 

Robust  Controller 

4.4815 

3.9349 

4 

Table  4  AO  Testbed  Disturbance  Rejection  Simulations. 


Clearly,  the  robust  controller  is  better  able  to  handle  disturbances  than  the 
classical  controller  is.  This  is  accomplished  in  fewer  iterations  as  well.  In  each  case,  the 
robust  controller  achieved  its  steady  state  wavefront  error  in  4  iterations,  whereas  the 
classical  controller  took  approximately  10  iterations  and  even  then  was  not  able  to 
completely  dampen  the  disturbance  input.  Similar  results  were  obtained  for  other  test 
frequencies.  The  additional  benefit  of  the  robust  controller  is  that  no  a  priori  knowledge 
of  the  disturbance  frequency  was  necessary.  It  handled  all  tested  frequencies  equally 
well.  It  is  this  ability  of  the  robust  controller  to  provide  robust  performance  in  the 
presence  of  uncertainty  that  makes  it  desirable  for  these  applications. 
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VII.  SUMMARY,  FUTURE  WORK,  AND  CONCLUSIONS 


A.  SUMMARY 

This  research  has  provided  an  important  contribution  to  the  field  of  control 
systems  for  adaptive  optics  systems.  It  has  demonstrated  the  design  of  a  robust  controller 
for  a  complex  AO  system.  While  not  the  first  use  of  Hw  control  for  AO  systems,  it  is  the 
first  analytical  work  for  a  space  based  AO  telescope  with  a  high  degree  of  complexity 
due  to  the  number  of  actuators,  sensors,  and  states.  To  accomplish  this,  an  innovative  use 
of  Zernike  polynomials  was  used  to  perform  model  reduction.  Zemike  polynomials  had 
previously  not  been  used  in  this  fashion  with  the  intent  of  facilitating  the  design  of  a 
robust  controller.  This  application  was  possible  due  to  the  fact  that  the  sensor  outputs  for 
the  AO  system  model  consisted  largely  of  Shack-Hartmann  wavefront  sensor  slopes. 
Zernike  polynomials  are  uniquely  suited  to  describe  this  type  of  data. 

An  analytical  model  that  originally  had  997  inputs  and  936  outputs  was  reduced 
to  21  inputs  and  21  outputs.  The  reduction  was  performed  by  a  Singular  Value 
Decomposition  approach  and  the  Zemike  polynomial  approach,  and  the  results  were 
compared.  Additionally,  a  Hankel  singular  value  reduction  was  performed  to  reduce  the 
number  of  states  in  the  system  from  332  to  240.  Both  input/output  reduction  approaches 
allowed  the  successful  synthesis  of  a  robust  controller.  The  controller  from  the  SVD 
reduced  model  had  a  closed  loop  control  bandwidth  between  20  and  24  Hz,  while  the 
controller  from  the  Zernike  reduced  model  had  a  closed  loop  control  bandwidth  between 
15  and  19  Hz.  Both  of  these  bandwidths  exceeded  the  required  10  Hz  bandwidth  that  the 
existing  classical  control  techniques  were  unable  to  achieve. 

While  classical  controllers  have  many  advantages  and  are  much  simpler  to  design 
and  implement,  the  ever  increasing  complexity  of  AO  systems  and  growing  number  of 
sensors  and  actuators  means  that  alternative  methods  need  to  be  developed.  The  work 
presented  here  has  demonstrated  that  robust  control  techniques  are  worthy  of  further 
research  for  implementation  in  these  types  of  systems. 
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Finally,  the  robust  controller  synthesis  design  techniques  used  for  the  complex 
analytical  model  were  experimentally  verified  on  an  AO  testbed.  While  the  AO  testbed 
system  was  much  simpler  than  the  analytical  space  telescope  model,  the  same  techniques 
and  algorithms  used  for  the  analytical  model  were  used  for  the  testbed. 

Both  controllers  achieved  similar  absolute  wavefront  correction  results;  however, 
the  robust  controller  was  able  to  achieve  the  steady  state  wavefront  error  in  only  4 
iterations,  whereas  the  classical  controller  took  approximately  25  iterations.  The  real 
advantage  of  the  robust  controller  was  seen  in  the  rejection  of  a  disturbance  input.  The 
robust  controller  was  essentially  impervious  to  the  disturbance  input  in  the  test  cases. 
This  was  possible  for  any  desired  sinusoidal  disturbance  frequency  without  any  a  priori 
knowledge  of  the  disturbance  frequency.  For  classical  control,  a  notch  filter  can  be  used 
for  disturbance  rejection.  However,  to  design  the  filter,  the  disturbance  frequency  has  to 
be  known  accurately.  The  superior  perfonnance  of  the  robust  controller  on  the  testbed 
versus  the  classical  controller  is  further  verification  that  the  techniques  implemented  in 
this  research  are  valid. 

B.  FUTURE  WORK 

A  number  of  issues  arose  during  this  research,  demonstrating  areas  where  further 
work  could  be  accomplished.  The  first  suggestion  for  additional  research  would  be  to 
experiment  with  different  numbers  of  Zernike  polynomials.  Twenty-one  polynomials 
were  used  in  this  research,  which  go  to  seventh  order  aberrations.  Using  more 
polynomials  could  improve  accuracy  of  the  wavefront  estimation.  The  tradeoff  would 
come  in  increased  complexity  and  possibly  longer  calculation  times  due  to  the  larger 
number  of  polynomials.  Additionally,  the  Zernike  polynomials  could  be  modified  for  an 
annular  ring  as  set  forth  in  [26]  to  account  for  the  central  hole  in  the  primary  mirror. 

Another  possibility  would  be  to  experiment  with  perfonning  a  state  reduction  at 
different  stages  in  the  process.  In  this  work,  a  state  reduction  was  performed  immediately 
after  the  inputs  and  outputs  were  reduced.  Alternative  approaches  would  be  to  perfonn 
state  reductions  on  the  original  full  plant,  the  augmented  2-input,  2-output  system  model, 
the  synthesized  robust  controller,  or  any  combination  of  multiple  state  reductions  at  these 
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stages.  An  important  consideration  for  the  segmented  space  telescope  model  is  that  it 
was  originally  given  in  modal  form  where  the  A  matrix  had  a  convenient  block  diagonal 
structure.  A  state  reduction  destroys  this  structure,  so  easy  identification  of  the  natural 
frequencies  is  no  longer  possible. 

Additionally,  the  robust  controller  synthesized  in  the  research  was  calculated 
based  on  a  direct  Hx  method.  Other  robust  controller  synthesis  techniques  exist  such  as 
the  /u  -synthesis  method,  which  might  yield  promising  results.  This  would  also  allow 
further  research  into  ways  of  introducing  uncertainty  into  the  segmented  space  telescope 
model. 

On  the  adaptive  optics  testbed,  refinement  of  the  wavefront  estimation  algorithms 
would  most  likely  provide  lower  absolute  error  measurements.  There  are  also  a  number 
of  other  classical  and  optimal  control  algorithms  that  could  be  implemented  and 
compared  on  the  testbed.  Upgrading  the  Basler  camera  on  the  testbed  used  for  wavefront 
sensing  to  a  higher  frames  per  second  rate  would  also  improve  the  system  performance. 

Eventually,  the  techniques  presented  here  could  be  implemented  on  a  testbed  with 
a  segmented  primary  optic.  The  Spacecraft  Research  and  Design  Center  has  such  a 
mirror  for  future  work.  Implementing  robust  control  techniques  on  this  mirror  could  be 
done  in  conjunction  with  the  research  currently  being  perfonned  at  the  SRDC  on  an 
improved  wavefront  sensing  technique  know  as  Redundant  Spacings  Calibration  [1]. 

C.  CONCLUSIONS 

This  research  has  contributed  to  the  state-of-the-art  in  the  following  areas: 

1.  An  Hx  robust  control  technique  has  been  used  for  surface  control  of  a 
flexible  space  telescope  model  including  structural  dynamics. 

2.  Robust  control  techniques  have  been  able  to  provide  a  desired  control 
bandwidth  that  was  not  possible  with  classical  control. 
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3.  A  new  model  reduction  technique  using  Zemike  polynomials  and  singular 
value  decomposition  has  been  developed.  This  technique  has  been  found  to 
provide  better  performance  to  reduce  the  wavefront  aberration  in  the  space 
telescope  model. 

4.  A  robust  control  technique  has  been  applied  experimentally  to  the  adaptive 
optics  testbed  of  the  Naval  Postgraduate  School  Spacecraft  Research  and 
Design  Center.  The  robust  control  technique  has  been  found  to  be  superior  to 
classical  control  in  response  time  and  disturbance  rejection. 

Hx  robust  control  techniques  offer  many  advantages  for  complex  AO  systems, 
particularly  those  that  are  subject  to  high  amounts  of  model  uncertainty,  such  as 
lightweight  segmented  space  telescopes.  As  future  system  become  more  complex,  and 
performance  requirements  become  more  stringent,  the  techniques  discussed  in  this 
dissertation  will  become  increasingly  important. 

It  is  important  to  note  that  just  because  a  robust  controller  can  be  designed  for  a 
particular  system,  this  does  not  necessarily  mean  that  it  should.  Classical  control 
methods,  optimal  control  techniques,  and  even  adaptive  controllers  all  offer  certain 
benefits.  The  ultimate  solution  will  most  likely  be  a  combination  of  approaches  to  reduce 
overall  system  risk,  guarantee  stability,  and  maximize  performance.  In  the  case  of  space 
systems,  the  extremely  high  development  and  launch  costs  require  multiple  levels  of 
redundancy  to  prevent  loss.  The  techniques  explored  in  this  research  will  undoubtedly  be 
applicable  to  future  designs. 
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